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Abstract 



QCD lattice simulations with 2+1 flavours (when two quark flavours 
are mass degenerate) typically start at rather large up-down and strange 
quark masses and extrapolate first the strange quark mass and then the 
up-down quark mass to its respective physical value. Here we discuss an 
alternative method of tuning the quark masses, in which the singlet quark 
mass is kept fixed. Using group theory the possible quark mass polynomials 
for a Taylor expansion about the flavour symmetric line are found, first for 
the general 1 + 1 + 1 flavour case and then for the 2 + 1 flavour case. 
This ensures that the kaon always has mass less than the physical kaon 
mass. This method of tuning quark masses then enables highly constrained 
polynomial fits to be used in the extrapolation of hadron masses to their 
physical values. Numerical results for the 2 + 1 flavour case conflrm the 
usefulness of this expansion and an extrapolation to the physical pion mass 
gives hadron mass values to within a few percent of their experimental 
values. Singlet quantities remain constant which allows the lattice spacing 
to be determined from hadron masses (without necessarily being at the 
physical point). Furthermore an extension of this programme to include 
partially quenched results is given. 
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1 Introduction 



The QCD interaction is flavour blind. Neglecting electromagnetic and weak in- 
teractions, the only difference between quark flavours comes from the quark mass 
matrix, which originates from the coupling to the Higgs field. We investigate here 
how flavour blindness constrains hadron masses after flavour SU{3) symmetry is 
broken by the mass difference between the strange and light quarks. The flavour 
structure illuminates the pattern of symmetry breaking in the hadron spectrum 
and helps us extrapolate 2 + 1 flavour lattice data to the physical point. (By 2 + 1 
we mean that the u and d quarks are mass degenerate.) 

We have our best theoretical understanding when all 3 quark flavours have 
the same masses (because we can use the full power of flavour SU{3) symme- 
try); nature presents us with just one instance of the theory, with mf/mf 25 
(where the superscript ^ denotes the renormalised mass). We are interested in 
interpolating between these two cases. We consider possible behaviours near the 
symmetric point, and find that flavour blindness is particularly helpful if we ap- 
proach the physical point, denoted by (mf*,mf*), along a path in the mf-mf 
plane starting at a point on the SU (3) flavour symmetric line (mf = = ttiq) 
and holding the sum of the quark masses = ^{m^ + + mf) = |(2mf + mf) 
constant pLj, at the value mg as sketched in Fig. [H The usual procedure (path) 



E 




^ 



Figure 1: Sketch of the path (red, solid line) in the mf-mf plane to the physical 
point denoted by {mf*,mf*). The dashed diagonal line is the S'f/(3)-symmetric line. 

is to estimate the physical strange quark mass and then try to keep it fixed, i.e. 
mf = constant, as the light quark mass is reduced towards its physical value. 
However on that path the problem is that the kaon mas^ is always larger than 
its physical value. Choosing instead a path such that the singlet quark mass is 

"'^In this article quark masses will be denoted by m, and hadron masses by M . 
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kept fixed has the advantage that we can vary both quark masses over a wide 
range, with the kaon mass always being hghter than its physical value along the 
entire trajectory. Starting from the symmetric point when masses are degenerate 
is particularly useful for strange quark physics as we can track the development 
of the strange quark mass. Also if we extend our measurements beyond the sym- 
metric point we can investigate a world with heavy up-down quarks and a lighter 
strange quark. 

The plan of this article is as follows. Before considering the 2-1-1 quark 
flavour case, we consider the more general 1 -|- 1 -f- 1 case in section O This 
also includes a discussion of the renormalisation of quark masses for non-chiral 
fermions. Keeping the singlet quark mass constant constrains the extrapolation 
and in particular it is shown in this section that flavour singlet quantities remain 
constant to leading order when extrapolating from a flavour symmetric point. 
This motivates investigating possible quark mass polynomials - we are able to 
classify them here to third order in the quark masses under the SU{3) and 5*3 
(flavour) groups. In section [3] we specialise to 2 -|- 1 flavours and give quark mass 
expansions to second order for the pseudoscalar and vector meson octets and 
baryon octet and decuplet. (The relation of this expansion to chiral perturbation 
theory is discussed later in section [5l) In section H] we extend the formalism to the 
partially quenched case (when the valence quarks of a hadron do not have to have 
the same mass as the sea quarks). This is potentially useful as the same expansion 
coefficients occur, which could allow a cheaper determination of them. We then 
turn to more specific lattice considerations in sections [6l [7] with emphasis on clover 
fermions (i.e. non-chiral fermions) used here. This is followed by section [HI which 
first gives numerical results for the constant singlet quark mass results used here. 
Flavour singlet quantities prove to be a good way of defining the scale and the 
consistency of some choices is discussed. We also investigate possible finite size 
effects. Finally in section [9] the numerical results for the hadron mass spectrum 
are presented in the form of a series of 'fan' plots where the various masses fan 
out from their common value at the symmetric point. Our conclusions are given 
in section (TDi Several Appendices provide some group theory background for this 
article, discuss the action used here and give tables of the hadron masses found. 

Mostly we restrict ourselves to the constant surface. However, in a few sec- 
tions, we also consider variations in (for example in the derivation of the quark 
mass expansion polynomials, section [231 the discussion of 0{a) improvement in 
section [231 and in section [13 where we generalise a constant fn^ formula). 

2 Theory for 1 + 1 + 1 flavours 

Our strategy is to start from a point with all three sea quark masses equal, 

K = ^d=^f = ^o^ (1) 



6 



and extrapolate towards the physical point, {m^* , m^* , mf*) , keeping the average 
sea quark mass 

m« = i(m:: + mj + mf), (2) 

constant at the value ttiq. For this trajectory to reach the physical point we 
have to start at a point where ttlq ^ ^mf*. As we approach the physical point, 
the u and d quarks become lighter, but the s quark becomes heavier. Pions are 
decreasing in mass, but K and r] increase in mass as we approach the physical 
point. 



2.1 Singlet and non-singlet renormalisation 

Before developing the theory, we first briefly comment on the renormalisation 
of the quark mass. While for chiral fermions the renormalised quark mass is 
directly proportional to the bare quark mass, = ZmiTiq, the problem, at least 
for Wilson-like fermions which have no chiral symmetry, is that singlet and non- 
singlet quark mass can renormalise differently [21 [3j§ 

= Z^{mq -m) + Z^Jn , q = u,d,s, (3) 

where niq are the bare quark masses, 

m = ^{m^ + rud + rus) , (4) 

is the non-singlet renormalisation constant, and Z^ is the singlet renormal- 
isation constant (both in scheme R). It is often convenient to re- write eq. 
as 

= Z^%mg + azm) , (5) 

where 

Z^ 

az = rm-l, ^m = (6) 



m 



represents the fractional difference between the renormalisation constants. (Nu- 
merically we will later see that this factor is ~ 0(1), and is thus non- negligible 
at our coupling.) This then gives 

m'' = z^%i + az)m . (7) 

This means that even for Wilson-type actions it does not matter whether we keep 
the bare or renormalised average sea quark mass constant. Obviously eq. (I7j) also 
holds for a reference point (mo,mo,mo) on the flavour symmetric line, i.e. 

= Z^mo = + (^z)mo . (8) 



^Perturbative computations showing this effect, which starts at the two-loop order, are given 
in S 
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Furthermore introducing the notation 

Snig = nig — , dniq = niq — fn, q = u,d,s , (9) 

for both renormahsed and bare quark masses, we find that 

Sm- = Z-^Smq. (10) 

So by keeping the singlet mass constant we avoid the need to use two different 
Zs and as we will be considering expansions about a flavour symmetric point, 
they will be similar using either the renormahsed or bare quark masses. (Of 
course the value of the expansion parameters will be different, but the structure 
of the expansion will be the same.) We shall discuss this point a little further in 
section 12.41 

So in the following we need not usually distinguish between bare and renor- 
mahsed quark masses. 

Note that it follows from the deflnition that 

Sniu + 6md + 6ms = , (11) 

so we could eliminate one of these symbols. However we shall keep all three 
symbols as we can then write some expressions in a more obviously symmetrical 
form. 



2.2 General strategy 

With this notation, the quark mass matrix is 



M 















-1 






V 





/ 






mSmu-Smd)\ 0-10 +i5mj -10 . (12) 



The mass matrix Ai has a singlet part (proportional to /) and an octet part, 
proportional to A3, As. We argue here that the theoretically cleanest way to 
approach the physical point is to keep the singlet part of Ai constant, and vary 
only the non-singlet parts. 

An important advantage of our strategy is that it strongly constrains the 
possible mass dependence of physical quantities, and so simplifles the extrapo- 
lation towards the physical point. Consider a flavour singlet quantity, which we 
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shall denote by X5, at a symmetric point (mo, mo, mo). Examples are the sea 1^ 
Xr = Tq^ , OT the plaquette P (this will soon be generalised to other singlet quan- 
tities). If we make small changes in the quark masses, symmetry requires that 
the derivatives at the symmetric point are equal 



dXs dX, 



dXs 
(9m, 



(13) 



So if we keep rriu + m^ + m^ constant, then any arbitrary small changes in the 
quark masses mean that Am^ + Am^ + Am^ = so 



AXc 



-/Xrriu + — — Am^ + — Am., 



drriu 



drrin 



0. 



(14) 



The effect of making the strange quark heavier exactly cancels the effect of making 
the light quarks lighter, so we know that Xs must be stationary at the symmet- 
rical point. This makes extrapolations towards the physical point much easier, 
especially since we find that in practice quadratic terms in the quark mass ex- 
pansion are very small. Any permutation of the quarks, such as an interchange 
M -H- s, or a cyclic permutation -u — )■ (i — )■ s — )■ m does not change the physics, 
it just renames the quarks. Any quantity unchanged by all permutations will be 
flat at the symmetric point, like X^. 

We can also construct permutation-symmetric combinations of hadrons. For 
orientation in Fig. |2] we give the octet multiplets for spin (pseudoscalar) and 



K^ds) ' 




K(us) 










I 






7i"(du)\ 


71° 


/7i^(ud) 


p-(du)\ 




K(su) 




K(sd) 


K{su) 





K(us) 




p"(ua) 



— *o_ 
K(sd) 



Figure 2: The octets for spin (pseudoscalar) and spin 1 (vector) mesons (plotted in 
the plane), r/g and (p^ are pure octet states, ignoring any mixing with the singlet 
mesons. 

spin 1 (vector) mesons and in Fig. [3] the lowest octet and decuplet multiplets for 

^There is no significance here to using rg or Tq ; however defining Xj, — Tq is more 
consistent with later definitions. 
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Y 




Figure 3: The lowest octet and decuplet for the spin ^ and for the spin | baryons. 

the spin ^ and for the spin | baryons (all plotted in the I3-Y plane). 

For example, for the decuplet, any permutation of the quark labels will leave 
the T,*^{uds) unchanged, so the is shown by a single black (square) point in 
Fig. m On the other hand, a permutation (such as u — )■ ci — )■ s) can change a 

00 A O O A 

o , ■ O O ■ O 

00 



Figure 4: The behaviour of the octet and decuplet under the permutation group 5*3. 
The colours denote sets of particles which are invariant under permutations of the quark 
flavours (red or filled triangles, blue or open diamonds and black or filled squares). 

A^^{uuu) into a A~{ddd) or (if repeated) into an Q~{sss), so these three particles 
form a set of baryons which is closed under quark permutations, and are all given 
the same colour red (triangle) in Fig. HI Finally the 6 baryons containing two 
quarks of one flavour, and one quark of a different flavour, form an invariant set, 
shown in blue (diamond) in Fig. HI 

If we sum the masses in any of these sets, we get a flavour symmetric quan- 
tity, which will obey the same argument we gave in eq. iHM for the quark 
mass (in) dependence of the scale tq. We therefore expect that the E*" mass 
must be flat at the symmetric point, and furthermore that the combinations 
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(Ma++ + Ma- + Mo-) and (Ma+ + Mao + M^,+ + M^*- + M=,o + Ms*-) will 
also be fiat. Technically these symmetrical combinations are in the Ai singlet 
representation of the permutation group 5*3. This is the symmetry group of an 
equilateral triangle, C^v This group has 3 irreducible representations, [6], two 
different singlets, Ai and A2 and a doublet E, with elements E~^ and E~. Some 
details of this group and its representations are given in Appendix |Xl while Ta- 
ble [1] gives a summary of the transformations. 





A, 


E 


A2 


Operation 




E+ E- 




Identity 
u ^ d 


+ 
+ 


+ + 
+ 


+ 


u ^ s 


+ 


mix 




d ^ s 


+ 


mix 




u ^ d s u 


+ 


mix 


+ 


u ^ s d u 


+ 


mix 


+ 



Table 1: A simplified table showing how the group operations of 5*3 act in the different 
representations: + refers to unchanged; — refers to states that are odd under the group 
operation. 



We list some of these invariant mass combinations in Table [2j The permuta- 



Pseudoscalar 




= i(M|+ + M|o + + M2_ + Mlo + M|.) 
= l«o + M„2j 


blue 


mesons 


x^ 


black 


Vector 


Xp- 


= 1{Mk*+ + Mk*o + Mp+ + Mp- + M^*o + Mk'-) 


blue 


mesons 


X<l>8 

X<p, 


= i(M^o + M^J 
= i(2M,o+M<^J 


black 


Octet 


Xn 


= i(Mp + M„ + M2+ + Ms- + Mho + M^-) 


blue 


baryons 


Xa 


= i(MA + Mso) 


black 


Decuplet 


Xa 


= i(MA++ +Ma~ +Mf,-) 


red 


baryons 




= \{Ma+ + Mao + Ms*+ + Ms*- + Mh.o + M^,-) 


blue 






= M2*o 


black 



Table 2: Permutation invariant mass combinations, see Fig. [H is a fictitious ss 
particle; ryg and 08 are pure octet mesons. The colours in the third column correspond 
to Fig. H 

tion group S3 yields a lot of useful relations, but cannot capture the entire struc- 
ture. For example, there is no way to make a connection between the A~^~^{uuu) 
and the A'^{uud) by permuting quarks. To go further, we need to classify phys- 
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ical quantities by SU{3) (containing the permutation group S3 as a subgroup), 
which we shall consider now. 

2.3 Taylor expansion 

We want to describe how physical quantities depend on the quark masses. To do 
this we will Taylor expand about a symmetric reference point 

{mu,md,ms) = {mo,mo,mo) . (15) 

Our results will be polynomials in the quark masses, we will express them in 
terms of m and 6mq of eq. (Q. The main idea is to classify all possible mass 
polynomials by their transformation properties under the permutation group 5*3 
and under the full flavour group SU{3), and classify hadronic observables in the 
same way. m and 6mq are a natural basis to choose as m is purely singlet and 
Snig is non-singlet. The alternative nig — niQ would be less useful as it contains a 
mixture of singlet and non-singlet quantities. 

The Taylor expansion of a given observable can only include the polynomials 
of the same symmetry as the observable. The Taylor expansions of hadronic 
quantities in the same SU{3) multiplet but in different S3 representations will 
have related expansion coefficients. (We will show examples of the latter, e.g. in 
eqs. (I3I1)-(I33]).) 

While we can always arrange polynomials to be in definite permutation group 
states, when we get to polynomials of 0{Smg) we find that a polynomial may be 
a mixture of several SU (3) representations, but the classification is still useful. In 
Table [3] we classify all the polynomials which could occur in a Taylor expansion 
about the symmetric point, eq. ( fT5l) . up to 0{6mg). 

Many of the polynomials in the table have factors of (m — mo). These polyno- 
mials drop out if we restrict ourselves to the surface of constant m = mo, leaving 
only the polynomials marked with a tick (/) in Table [31 At 0{m^) there are 
k + 1 independent polynomials needed to describe functions on the constant m 
surface (the polynomials with the ticks), but ^{k + l){k + 2) polynomials needed 
if the constraint m = constant is dropped (all polynomials, with and without 
ticks). Thus the advantage of working in the constant m surface increases as we 
proceed to higher order in m^. 

Since we are keeping m constant, we are only changing the octet part of the 
mass matrix in eq. (fT2l) . Therefore, to first order in the mass change, only octet 
quantities can be affected. SU{3) singlets have no linear dependence on the quark 
mass, as we have already seen by the symmetry argument eq. f|T^ . but we now 
see that all quantities in SU{3) multiplets higher than the octet cannot have 
linear terms. This provides a constraint on the hadron masses within a multiplet 
and leads (as we shall see) to the Gell-Mann-Okubo mass relations [71 [S] . 

When we proceed to quadratic polynomials we can construct polynomials 
which transform like mixtures of the 1, 8 and 27 multiplets of SU{3), see Table |3l 
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Polynomial 




5-3 




SU{3) 






1 


/ 


Ai 


1 








(m — mn) 




Ai 


1 










/ 


E+ 




] 








/ 


E- 




] 






(m — TTln)^ 

V U / 




A^ 


1 








(m — ■mQ)6ms 




E+ 




i 






(fn — Tnr\)(5m,, — 8171,1) 




E- 










6m'^ + 6m'd + Smj 


/ 


Ai 


1 




27 




Sdml — {drriu — Snid)'^ 


/ 


E+ 




i 


27 






/ 


E' 




] 


27 




im — tuqY 




A 


1 
i 












F+ 


c 


2 
5 






{ qm r^\'^ ( Sim Sim t\ 

yin irLQ) yuiiLu 'J' 'I'd) 
(m — mo)(5m^ + + (5m^) 




11/ 


i 

c 


2 










1 




27 




(m — mo) [35m^ — (5m„ — SmdY] 




E+ 




] 


27 




(m — mo)5ms((5mrf — 5m„) 




E- 




] 


27 




6mu6mdSms 


/ 


Ax 


1 




27 


64 


6ms{Sml^ + 5rn?^ + 5m1) 


/ 


E+ 




] 


27 


64 


{drriu — 6md){Sml^ + 5rn?^ + 5m1) 


/ 


E- 




] 


27 


64 


{6ms - 6mu){Sms - 6md){Smu - Snid) 


/ 


A2 




10 To 




64 



Table 3: All the quark- mass polynomials up to 0{mg), classified by symmetry prop- 
erties. A tick (/) marks the polynomials relevant on a constant m surface. These 
polynomials are plotted in Fig. [6l If we want to make an expansion valid when fn 
varies, then all the polynomials in the table (with and without ticks) are needed. 



Further representations, namely the 10, 10 and 64, first occur when we look at 
cubic polynomials in the quark masses, see again Table El 

In a little more detail, constructing polynomials with a definite 5*3 classifica- 
tion is fairly straightforward, we have to see what happens to each polynomial 
under simple interchanges (e.g. u d) and cyclic permutations (e.g. m — )■ s, 
s — )■ d — )■ u). The ^3 column of Table [3] is easy to check by hand. The 
SU{3) assignment of polynomials is less straightforward. Only the simplest poly- 
nomials belong purely to a single SU{3) multiplet; most polynomials contain 
mixtures of several multiplets. The non-singlet mass is an octet of SU{3), so 
quadratic polynomials in 6mq can contain representations which occur in 8 (g) 8, 
cubic polynomials representations which occur in 8 08® 8. We can find out which 
representations are present in a given polynomial by using the Casimir operators 
of SU{3) P [in]. That operator was programmed in Mathematica, and used to 
analyse our polynomial basis. Some more details are presented in Appendix [B] 
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(in section [R2|) . The results of the calculation are recorded in the SU{3) section 
of Table [31 

The allowed quark mass region on the fn = constant surface is an equilateral 
triangle, as shown in Fig. O Plotting the polynomials of Table E] across this 
triangular region then gives the plots in Fig. El where the colour coding indicates 
whether the polynomial is positive (red) or negative (blue). 

As a first example of the use of these tables, consider the Taylor expansion 
for the scale ro/a up to cubic order in the quark masses. As discussed previously, 
this is a gluonic quantity, blind to flavour, so it has symmetry Ai under the S3 
permutation group. Therefore its Taylor expansion only contains polynomials of 
symmetry Ai. If we keep m, the average quark mass, fixed, the expansion of ro/a 
must take the form 

— = a + I3i5m^ + 5rn?^ + 5rn?) + 7 6mu6mdSms , (16) 
a 

with just 3 coefficients. Interestingly, we could find all 3 coefficients from 2 + 1 
data, so we would be able to predict 1 + 1 + 1 flavour results from flts to 2 + 1 data. 
This is common. If we allow m to vary too, we would need 7 coefficients to give a 
cubic fit for tq (all the Ai polynomials in Table [31 both ticked and unticked). This 
point is further discussed in section 14.41 If we did not have any information on 
the flavour symmetry of tq we would need all the polynomials in Table [3l which 
would require 20 coefficients. 



2.4 0(a) improvement of quark masses 

Before classifying the hadron mass matrix, we pause and consider the 0{a) im- 
provement of quark masses. (If we are considering chiral fermions, we have 'auto- 
matic 0{a) improvement', see e.g. [11] for a discussion.) In writing down expres- 
sions for bare and improved quark masses, it is natural to expand about the chiral 
point, all three quarks massless, which means setting mo = in the expressions 
in Table [3l Later, when we consider lattice results, we want to expand around a 
point where we can run simulations, so we will normally have a non-zero rriQ. 

Improving the quark masses requires us to add improvement terms of the 
type am^ to the bare mass. We can add S'?7(3)-singlet improvement terms to 
the singlet quark mass, SU (3)-octet improvement terms to the non-singlet quark 
mass. We are led to the following expressions for the improved and renormalised 
quark masses 

m"* = Z^[m + a [bifrf + b2{Sml + 6ml + 6ml)]] 
6mf = Z^^ [6ms + a {b3m6ms + h{36ml - {6m^ - 6md)^)}] , (17) 

together with = Z^Vm, eqs. (El), ([6]). We have improved by adding the 
two possible singlet terms from the quadratic section of Table [31 and improved 
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Figure 5: The allowed quark mass region on the m = constant surface is an equilateral 
triangle. The black point at the centre is the symmetric point, the red star is the 
physical point. 2 + 1 simulations lie on the vertical symmetry axis. The physical point 
is slightly off the 2 + 1 axis because > ruu- 




Figure 6: Contour plots of the polynomials relevant for the constant m Taylor ex- 
pansion, see Table El A red(dish) colour denotes a positive number while a blue(ish) 
colour indicates a negative number. If rriu = nid (the 2 + 1 case), only the polynomials 
in the Ai and columns contribute. Each triangle in this figure uses the coordinate 
system explained in Fig. [5l 
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6mf by adding the two possible octet polynomials. Note that if we keep fn 
constant, we only need to consider the improvement terms 62 and 64. The bi and 
63 terms could be absorbed into the Z factors. Simplifications of this sort are 
very common if m is kept fixed. 

We get expressions for the u and d quark mass improvement by flavour- 
permuting eq. ( !T7|) 

(5m^ = [5mu + a{b3m6mu + h{36ml- {5ms - drndY)}] 

6m^-6m^ = Z^[5mu-5md (18) 

+a {h^mi^dmu — Smd) + 6645ms (^m^ — 5m„)}] . 

The improvement terms for SrUu — Srud are proportional to the two E~ , SU{3)- 
octet, quadratic polynomials. (We have to use the identity, eq. f|TT]) . to bring the 
result to the desired form - which will often be the case in what follows.) 

Table |3] is based purely on flavour arguments, we would hope that all the 
results are true whether we use bare or renormalised quantities, and also in- 
dependently of whether we work with a naive bare mass, or a bare mass with 
0{a) improvement terms. Let us check if this is true. The first thing we need 
to know is whether the zero-sum identity eq. (fTTj) survives renormalisation and 
improvement. Using the previous equations we find 

5m^ + Sm^ + Smf = Z^^ [{Srriu + Sma + Sm^) 

+a {bsfn^drriu + Snid + 6ms) 
+64(5m„ + 6md + 6msY}] 
= 0, (19) 

showing that eq. (fTTj) is not violated by improvement or renormalisation. 

The next point we want to check is if the symmetry of a polynomial depends on 
whether we expand in terms of improved or unimproved masses. As an example, 
let us look at the quadratic polynomial 

Sm^iSm^ - 5m«) , (20) 

which has permutation symmetry , and SU{3) content octet and 27-plet. 
Expanding to first order in the lattice spacing a we find 

(5mf (5m^ - 5m^) = (Z;^^)^ [5ms((5mrf - 5m„) 

+a {2bsfnSms{Smd — Smu) (21) 
-I- 2bA{5mu - 5md){Sml + Sm^ + 6ml)}] . 

The mass improvement terms have generated two extra cubic polynomials, but 
they are both polynomials with the same symmetry as the initial polynomial. 
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The same holds for the other quadratic terms. This shows that Table |3] applies 
both to improved and unimproved masses. 

Thus our conclusion is that the flavour expansion results are true whether we 
use bare or renormalised quantities, and also independently of whether we work 
with a naive bare mass, or a bare mass with 0(a) improvement terms. 

Finally we compare these results with those obtained in [12], to see whether we 
can match the 4 improvement terms found in eq. f|T71) to the 4 terms introduced 
there, namely 



rriq + {rm - l)m (22) 



+a <^ bmml + Sbmiriam + irmdm - bm)m^ + "^irmdm - bm)m^ 



where w? = |(m^ + m^ + m^). At first this looks different from eq. ( IT7|) . but this 
is just due to a different choice of basis polynomials. The quadratic polynomials 
in eq. f l22p are simple linear combinations of those in eq. f ll7p . 
From eq. (fT7|) we have 

mf = 6m^ + m'' (23) 
= [nis + — ^)rn + a [b^rnSms + b^^SSml — (^m^ — SnidY) 

+rmbirrf + rj32{5ml + 6ml + ^'^d)}] ^ 
so we now equate the terms to those in eq. (1221) . We must first re- write 



so 



3Sml — (Sniu — SnidY = 6 



,2 o™, 



2mm= — + 2m^ 



(24) 



+a < b-^m{ms — m) + 664(7/1^ — 2mms — m? + 2m ) 







+rmbim'^ 


+ 3rmb2{m'^ 








which gives the results 












bm 


= 664 




bi 


= 3dm 


+ 


dm 


bm 
dm 


= 

= 362 


-464 


62 

or , 

03 


— 3bm 


+ 


2bm 


dm 


= 


-62 


64 


- ifc 

g 







(25) 



(26) 



2.5 SU(3) and S3 classification of hadron mass matrices 

In eq. ( fT2l) we split the quark mass matrix into a singlet part and two octet parts. 
We want to make a similar decomposition of the hadron mass matrices. We start 
with the decuplet mass matrix because it is simpler than the octet mass matrix. 
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2.5.1 The decuplet mass matrix 



The decuplet mass matrix is a 10 x 10 diagonal matrix. From SU{3) group algebra 
we know 

10®I0 = 1©8©27©64. (27) 

The singlet matrix is the identity matrix, the octet representation contains 2 
diagonal matrices (A3 and Ag), the 27-plet has 3 diagonal matrices, and the 64- 
plet includes 4 diagonal matrices, see Fig. [71 This gives us a basis of 10 diagonal 



Figure 7: An illustration of the octet, 27-pIet and 64-plet representations of SU{3). 
The number of spots in the central location gives the number of flavour-conserving 
operators in each multiplet. In the octet, the 2 operators form an E doublet of the 
permutation group. In the 27-plet the 3 operators are an Ai singlet and an E doublet. 
In the 64-plet the centre operators are an Ai singlet, an E doublet and an A2 singlet. 

matrices, into which we can decompose the decuplet mass matrix. 

We can use the Casimir operator to project out the diagonal matrices in a 
particular SU{3) representation (see Appendix IB. 31 for a fuller discussion). As 
an example of a matrix with pure octet symmetry, we can take the operator 2/3. 
(We have multiplied by 2 simply to avoid having fractions in the matrix.) Since 
we know the isospins of all the decuplet baryons, we can write down 

A" A° A+ A++S* 

/ -3 
0-100 

10 

3 











\ 







+ 















































































-2 









































2 




















-1 




















1 




















) 



(28) 



-3-113 
-2 2 
-1 1 
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where we have used a more compact notation to record the diagonal elements 
on the right-hand side. The entry in the A~ column of the matrix is —3, so on 
the right-hand side we put a —3 in the position of the A~ in the usual decuplet 
diagram, and so on. By considering the reflection and rotation symmetries of the 
right-hand side of eq. fj28|) we can see that this matrix corresponds to the basis 
element E~ of the doublet representation of 5*3. 

In Fig. [8] we show all 10 diagonal matrices, in this compact notation. These 

Ai E" A2 

1111 
1 111 
1 1 
1 

8 



3-1-1 3 
27 -1 -3 -1 
-1 -1 
3 



1 -1 -1 -1 
000 
1 1 
2 



-377-3 
-5 -5 
-2 -2 

6 



-3-1 1 3 
-2 2 
-1 1 




-3-11 3 
3 0-3 
4 -4 




64 



2-3-3 2 
-3 12 -3 
-3 -3 
2 



-100-1 

3 3 
-3 -3 
2 



-12-2 1 
1 -1 
-1 1 




0-110 
1 -1 

-1 1 





Figure 8: The matrices for projecting out decuplet mass contributions of known 
symmetry - see eq. ()28p for an explanation of the notation. 

matrices are orthogonal, in the sense 

Tr[r,r,]=0 if a^b, (29) 

(where is any of the matrices of Fig. |8]) so they can be used to project out 
mass combinations which have simple quark mass dependencies, see Fig. IH and 
Table H 

Let us now give some examples of mass formulae. First we look at the singlet 
of the decuplet mass matrix. Because we are keeping m = constant only the 
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Table 4: Decuplet mass matrix contributions, classified by permutation and SU{?>) 
symmetry, see Fig. [8j 

terms with ticks in Table |3] contribute. This gives from Table IH 

Ma- + Mao + Ma+ + Ma++ 
+Ms*- + Ms.o + Ms.+ + Mh.- + Mh.o + Mq- 
= lOMo + Bx {6ml + 6ml + 6ml) + Ci 6mjmd6ms . (30) 

This equation being a singlet has the same form as for rg/a, eq. f|T6|) . 

As a further example for the 27-plet component of the decuplet mass matrix, 
we see from Table H] that there are three mass combinations which transform as 
27-plets, giving three related mass relations 

3Ma- - Mao - Ma+ + 3Ma++ 

-Ms*- - 3Me,o - Ms*+ - Mh*- - Mh.o + 3Mn- 

= 627 [6'ml + 6ml + 6ml] + 9c276mu6md6ms (31) 
-3Ma- + 7Mao + 7Ma+ - 3Ma++ 
-5Ms.- - 5Ms.+ - 2Mh*- - 2Mh*o + 6Mf^- 

= 627 [3(5mg - {6mu - ^m^)^] + 3c276ms {6ml + 6ml + 6ml) (32) 
-3Ma- - Mao + Ma+ + 3Ma++ 
+3Ms*- - 3Ms.+ + 4M=,- - 4M=,o 

= 2h27{6md - 6mu)6ms + C27{6mu - 6md) {6ml + 6ml + 6ml) . (33) 

The coefficients in eqs. fl3T]) - fl33|) are connected, they all involve just one quadratic 
parameter, 627, and one cubic parameter, C27. We now want to explain the differ- 
ent numerical coefficients in front of these parameters. These can be checked by 
considering some simple symmetry limits. First consider the isospin limit, equal 
masses for the u and d quarks, (5m„ — )■ 6mi, 6md — )■ 6mi, 6ms — )■ —26mi (from 
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eq. (mD). In this limit, eq. §^ reduces to = 0, while eqs. §^ both 

become 

AM A - 5Me* - 2Mh. + 3Mn = Qb276mf - IScar^mf . (34) 

To include eq. fl33|) in our checks, we can take the ?7-spin limit, rus — )■ nid, i.e. 
6ms 6md, Srriu —26md. In this limit all decuplet baryons with the same 
electric charge would have equal mass, because they would be in the same f/-spin 
multiplet, so Mq- M^-^M^*- Ma-,Ms*- Ma- and similarly for the 
other charges. Now, all three equations become identical, 

4Ma- - 5Mao - 2Ma+ + 3Ma++ = 6b276ml - IScay^mj^ , (35) 

which again confirms that the numerical coefficients in eqs. (13T])-( 133|) are correct. 

Finally note that we can find all the coefficients in these equations from a 2 + 1 
simulation, and use them to (fully) predict the results of a 1 + 1 + 1 simulation. 

2.5.2 The octet mass matrix 

We can analyse the possible terms in the octet mass matrix in the same way 
as we did for the decuplet. We first consider the baryon octet. Using the same 
technique as for the decuplet mass matrix we find the results given in Table |5l 
However there is a complication in the octet case which we do not have in the 
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Table 5: Mass matrix contributions for octet baryons, classified by permutation and 
SU{3) symmetry. Note that the first two octet quantities (the 8a) are proportional to 
the hypercharge Y and to isospin I3, respectively. 

decuplet, caused by the fact that we have two particles (the A and S°) with the 
same Y and I3 quantum numbers. If 7^ these states mix. There are 
interesting connections between the elements of the A/S*^ mixing matrix and the 
splittings of the other baryons, but since in this article we are concerned with 
2 + 1 simulations, where this mixing does not arise, we will not discuss this further 
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here. We can however pick out several useful mass relations which are unaffected 
by A/S° mixing 

Mn + Mp + Ma + Ms- + Mso + Ms+ + M^- + M^o 

Mn + Mp- 3Ma + Ms- - SMso + + M^- + M^o 
= b27{6ml + 6ml + 6ml). (36) 

At order we meet some quantities in the baryon octet masses (the 10 and 
10 combinations) which can not be deduced from 2 + 1 flavour measurements - 
though valence 1 + 1 + 1 on a2 + l background would be a possible method of 
estimating these quantities. 

One early prediction concerning hyperon masses was the Coleman- Glashow 
relation [13] 

Mn-Mp- Ms- + Ms+ + Ms- - Mho ^ . (37) 

Deviations from this relation are barely detectable, using a recent precision mea- 
surement of the S° mass [I4] gives the value [15], M„ — Mp — Ms- + Ms+ + 
Ms- — Mso = —0.29 ± 0.26 MeV. The original Coleman-Glashow argument 
showed why the leading electromagnetic contribution to this quantity vanishes 
(in modern terms, the leading electromagnetic mass contributions are unchanged 
by the operation -H- s because the s and d quarks have the same charge, but the 
quantity in eq. (1371) is odd under this operation). To understand the smallness 
of the Coleman-Glashow quantity we also need to explain why the contribution 
from flavour SU{3) breaking due to quark mass differences is small. The mass 
combination appears in Table [5] as one of the A2 quantities. We can understand 
the success of the Coleman-Glashow relation by noting that the only polynomial 
in Table m with A2 symmetry is 0{6mg), so that the predicted violation of the 
Coleman-Glashow relation is 

Mn- Mp- Ms- + Ms+ + Me- - M^o 
= Cio{6i7is - 6mu){6ms - 6i7id){6mu - 6md) . (38) 

The polynomial is zero if any pair of quarks have the same mass, so we would 
need to measure the masses of baryons in a 1 + 1 + 1 setting to determine cio and 
predict the violation of the Coleman-Glashow relation. 

Turning now to the mesons, both the pseudoscalar and vector meson octet 
have a similar mass matrix, so they do not have to be considered separately. In 
Table |6] we give the mass matrix contributions for the octet mesons, classified by 
permutation and SU (3) symmetry. 

Some contributions allowed for baryons. Table El are absent for mesons be- 
cause they would violate charge conjugation, giving (for example) different masses 
to the and K~. In particular, there are two octets, 8a and 8;,, in the baryon 
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Table 6: Mass matrix contributions for octet mesons, classified by permutation and 
SU{3) symmetry. 

table, Table O but only the 8b octet is permitted in Table |6l When we write 
down the mass formulae, this will mean that the baryon mass formula will have 
two independent terms linear in the quark mass, but the meson mass formula 
will only have a single linear term. 

3 Theory for 2 + 1 flavours 

If we take any mass relation from the previous 1 + 1 + 1 section, and put m„ = 
ma = nil we will get a valid mass relation for the 2 + 1 case. In the 2 + 1 case 
only one variable is needed to parametrise the symmetry breaking, since from 
eq. (dH), 

Snig = —2Smi , (39) 

where 

5mi = mi — m . (40) 

In the 2 + 1 case we can simplify the mass matrix Tables HI O [61 The E~ and 
A2 matrices are not needed, their coefficients are always proportional to m^ — md, 
which we are now setting to zero. In the higher representations (27-plet and 64- 
plet) only one linear combination of the Ai and E~^ matrices appears in the 2 + 1 
case (it is the linear combination which does not split particles within an isospin 
multiplet). Therefore, in this section we just need the simplified matrix Tables [TJ 
[HlandEl 

In the 2 + 1 limit the decuplet baryons have 4 different masses (for the A, 
S*, S*, and Q). Similarly, for the octet baryons there are also 4 distinct masses, 
(iV, A, Z1,S); and for octet mesons, 3 masses. In the meson octet the K and 
K must have the same mass, but there is no reason why the N and S (which 
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Table 7: Decuplet mass matrix contributions for the 2 + 1 case, classified by SU (3) 
symmetry. Compare with Table [H 
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Table 8: Mass matrix contributions for octet baryons for the 2 + 1 case, classified by 
SU{3) symmetry. Compare with Table [5l 

occupy the corresponding places in the baryon octet, see Figs. [2] and [3]), should 
have equal masses once flavour SU{3) is broken. 

Again we have the singlet quantities Xs which are stationary at the symmetry 
point as given in Tabled but which now simplify to give Table fTOl In the notation 
we have now assumed isospin invariance, so that for example Ma = M^++ = 
Ma+ = Mao = Ma- . (The corresponding mass values we use in this article are 
given in section |9l) 

Furthermore this can obviously be generalised. Let us first define the quark 
mass combinations = {mi + 2ms)/ 3 and m^ = {mi + m^j2. Then mi + m.^ = 
2m and mi + 2mK = 3m are constants on our trajectory and so 6mi + (5m^ = 
and also 6mi + 26mK = 0. For example, this means that any functions of the 
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Table 9: Mass matrix contributions for octet mesons for the 2 + 1 case, classified by 
SU{3) symmetry. Compare with Table [6l 
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Pseudoscalar 
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x„ 
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Table 10: Permutation invariant mass combinations, see Fig. [H i;^<j is a fictitious ss 
particle; r/g a pure octet meson. 



form 

2f{mK) + /(m/) or g{ms) + 25((mi) or /i(m^) + h{mi) , (41) 

will also have zero derivative at the symmetric point. These equations generalise 
the meson sector of Table [IDl 

In Table [11] we present the 2 + 1 baryon decuplet results corresponding to 
Table m or eqs. fl30|) - fl33|) and the baryon octet from Table [5] or eq. fl36|) . Similarly 
the mesons are given in Table [121 Particular combinations chosen to remove the 
unknown M^g , M^^ masses are given in Table [131 

We can see how well this works in practice by looking, for example, at the 
physical masses of the decuplet baryons. If we consider the physical values of the 
four decuplet mass combinations in Table [Til and using mass values given later 
in Table [TTJ we get 

4Ma + 3Ms. + 2Ms* + Mq = 13.82 GeV singlet oc {5mif 

-2Ma + Mh. +Mn = 0.742 GeV octet oc 5mi (42) 

4Ma - 5Ms* - 2Mh* + 3Mq = -0.044 GeV 27 - plet oc 5mf 

-Ma + 3Ms. - 3Mh. + = -0.006 GeV 64 - plet oc 5m] , 

with a strong hierarchy in values, corresponding to the leading term in the Taylor 
expansion. Each additional factor of 5mi reduces the value by about an order 
of magnitude, the 64-plet combination is more than 2000 times smaller than the 
singlet combination. This suggests that the Taylor expansion converges well all 
the way from the symmetric point to the physical point. (Though of course it is 
possible that the singlet and octet curvature terms are larger than those in the 
27 and 64.) Unfortunately it may be very difficult to see a signal in the 64-plet 
channel. 
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Table 11: Baryon mass combinations classified by SU{3) representation, in the 2 + 1 
case. 
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Table 12: Meson mass combinations classified by SU{3) representation, in the 2 + 1 
case. Octet-singlet mixing is not taken into account. 



SU{3) Mass Combination Expansion 
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Table 13: Meson mass combinations free from mixing problems, classified by SU{3) 
representation. These combinations have been chosen to eliminate the rjs and 4>s states, 
so they now contain mixtures of different SU{3) representations. 
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We can now 'invert' the results in Tables [TT], [12] to give the expansion for 
each hadron mas^ from the symmetry point m = mo- This inversion is made 
easier by orthogonality relations between the different SU{3) representations. We 
can simply read off the answers from the Tables U\ M and M This leads to the 
constrained fit formulae 

= + 2a6mi + {f3o + 2f3i)6m^ 
Ml = -a6mi + {f3o + 5f3i + 9f32)6m^ 

Ml = Ml -2a5mi + {l3o + Q(5i + 12(52 + I3^)6m^i, (43) 



Mp = Mo + 2a5mi + {(3o + 2(3i)5mf 
Mk* = Mo - a5mi + (/3o + 5/3i + 9l32)Smf 
M^3 = Mo -2aSmi + i(3o + 6(3, + 12(32 + (3s)5m^, (44) 



Mjv = Mo + 3AiSmi + (Bq + 3Bi)5mf 

Ma = Mo + 3A26mi + {Bq + 6Bi - 3^2 + 9Bi)Sm^ 

Ms = Mo-3A25mi + {Bo + 6Bi + 3B2 + 9Bs)5mf 

Mh = Mo-3{Ai- A2)6mi + {Bo + 9Bi-3B2 + 9B3)6mf , (45) 



Ma = Mo + 3A5mi + {Bq + 3Bi)5m1 

Me. = Mo + + (5o + 65i + 9B2)Smf 

Ms* = Mo-3A5mi + {Bo + 9Bi + 9B2)5mf 

= Mo - 6A5mi + (Bo + 12Bi)dmf . (46) 

(The values of the constants are obviously different for each octet or decuplet.) 
We see that the linear terms are highly constrained. The decuplet baryons have 
only one slope parameter, A, while the octet baryons have two slope parameters, 
Ai, A2. Mesons have fewer slope parameters than octet baryons because of 
constraints due to charge conjugation, leaving again just one slope parameter. 

The quadratic terms are much less constrained; indeed only for the baryon 
decuplet is there any constraint. The coefficients of the 6mf terms appear com- 
plicated; there seem to be too many for the meson and baryon octets, eqs. (H3ll 
- f l45|) . In the next section, we generalise these formulae to the case of differ- 
ent valence quark masses to sea quark masses or 'partial quenching' when this 
choice of coefficients becomes relevant. Note that not all the coefficients can thus 
be determined: for rjs and (j)s the /Ss-coefficient cannot be found from partially 
quenched results. 

Alternatively, of course, the niethods of Appendix IB .41 can be directly applied. 
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An ss meson state has charge 0, isospin 0, the same quantum numbers as 
a {uu + dd) meson. In the real world we should therefore expect that / = 
(isoscalar) mesons will always be mixtures of ss and U, to a greater or lesser 
extent. The mixing has been investigated in detail in [16]. 

On the lattice we can remove the mixing by just dropping disconnected con- 
tributions, and only keeping the connected part of the meson propagator. (In fact 
this calculation is easier and cheaper than the full calculation.) Theoretically the 
resulting pure ss meson is best treated in the context of partial quenching. We 
can get its mass formula simply by making the substitution for the valence quark 
mass 6fis —25/1; in the mass formula for the partially quenched 1 = 1 mesons 
as described in section HJ eq. ( [57|) to give 



(These are, as previously mentioned, the masses that can be (easily) found from 
lattice simulations rather than M^g and M^^.) The pseudoscalar ss state, the rjs 
or 'strange pion', does not correspond to any real-world state, but the vector ss 
meson, the 0^, is very close to the real- world 0. Phenomeno logically, the observed 
fact that the almost always decays to KK (84 %) rather than pvr [17j is best 
explained by saying that the is almost purely ss. 

As eqs. fH3|) - fH6|) have been derived using only group theoretic arguments, 
they will be valid for results derived on any lattice volume (though the coefficients 
are still functions of the volume). 

Finally there is the practical question of whether fits should be against the 
(light) quark mass or alternatively against the pseudoscalar pion mass. In Ap- 
pendix O we argue that 'internally' at least the fits should be made against the 
quark mass. Of course this is only a useful observation when quadratic or higher 
terms are involved. To leading order eqs. (H^ - (l46|l can be written as 



(together with = — Mq = cmSttii in the case of pseudoscalar mesons, 
eq. fH3|) ). The coefficients cm can be found from these equations. Thus an 
expansion in 6mi is equivalent to an expansion in 6M or 5M^. 

4 Partial quenching 

In partial quenching (PQ) measurements are made with the mass of the valence 
quarks different from the sea quarks. In this case the sea quark masses mi, rris 
remain constrained by m = constant, but the valence quark masses /i;, fis are 
unconstrained. We define 




Mq - Aa5mi + (/3o + 8/3i)5mf . 



(47) 



6M = M - Mo = cmSttii , 




= jjLq — m , q = l,s . 



(49) 
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When fi ^ m (i.e. return to the 'unitary hne') then the following results collapse 
to the previous results of eqs. (H3|) - P6l) . Here we sketch some results, see [18] 
for more details. 

In Table [3] we have often used the identity drriu + 6md + Snis = to simplify 
the symmetric polynomials. Since we are not going to keep /iM + /id + /is constant, 
we write out our basis polynomials in full, without the constraint ^ 6 fig = 0, in 
Table [HI We can check that the polynomials in Table [H] reduce to multiples of 
those in Table |3] when the identity ^ Sfig = is applied. For example, 26 fig — 
Sfj,u — Sfid = 'iSfis — {6fiu+6fid+Sfj,s), so when the zero-sum identity is applied, and 
fiq — )• niq the polynomial 26fis—SfJ'u—SfJ'd in Table [H] reduces to the corresponding 



Polynomial 


S3 


SU{3) 


1 


Ai 


1 


6 flu + 6 fid + Sfis 


Ai 


1 


26fis - 6 flu - 6 fid 
5 flu - 6 fid 


E+ 
E- 


8 
8 


{6 flu + Sfid + SfisY 
{6 flu + 5 fid + 5fis){25fis - 5 flu - 5 fid) 
{6 flu + Sfid + 5fis){5fiu - 5 fid) 


Ai 
E+ 
E- 


1 

8 
8 


(Sfis - Sfiu)"^ + {SfXs - SfXd)"^ + {5fiu - 5fid)^ 
{5fis - 5fiu)^ + {Sfis - Sfid)^ - 2{5fiu - 5fid)^ 
{Sfis - Sfiu)"^ - {Sfis - Sfid)"^ 


Ai 

E+ 
E- 


1 27 
8 27 
8 27 


6ml + 6mj + 6m'j 


Ai 


1 27 



Table 14: All the quark-mass polynomials needed for partially quenched masses, 
classified by symmetry properties. The table includes entries up to 0{fig). 

polynomial 6ms in Table |3l 

4.1 PQ decuplet baryons 

In the partially quenched case we know that the hadron mass formulae should 
have an SU{3) symmetry for interchanging the sea quarks, and another SU{3) 
symmetry for operations on the valence quarks. The sea quark symmetry will 
always be singlet, the valence quark terms for a hadron in the octet can be in any 
representation which occurs in 8® 8 and the decuplet can be in any representation 
which occurs in 10 ® 10. 

Let us see what sort of mass relations symmetry allows us, taking the decuplet 
baryons as our example. Starting with linear terms in the quark masses, we can 
form two polynomials of the valence masses, a singlet combination {26fii+6fis) and 
an octet combination with E'^ symmetry, {6 fig — 6fii). These are the polynomials 
on line 2 and line 3 of Table [TH (A first-order term in the sea quark masses 
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is ruled out because we are keeping 2mi + constant.) We can read off the 
coefficients each polynomial must have by looking for the Ai singlet and 
octet entries in Fig. [HI or in Tables H] or [7l Thus a singlet polynomial must have 
the same flavour coefficients for every baryon and the octet polynomial must have 
a coefficient proportional to the hypercharge. 
So, at first sight we would expect 

Ma = Mo + ai{26fii + 6fis) - a2{6fis - Sfii) 

Ms* = Mo + ai{2Sfxi + SfXs) 

Mh* = Mo + ai{26jj,i + 6jj,s) + a2{6jj,s - 6jj,i) 

Mo = Mo + ai{26fii + 6fis) + 2a2{6fis-Sfii), (50) 

with no connection between ai and 02 • However, it is clear that the A mass 
cannot know anything about the strange valence mass, and the Q mass must 
similarly be independent of These constraints are both satisfied if 

ai = a2 = A , (51) 

giving us a leading-order formula 

Ma = Mo + 3A6fii 

Ms* = Mo + A{26fii + 6ns) 

Ms* = Mo + A{6fii + 26fis) 

Mn = Mo + 3A6fis. (52) 

We can continue this procedure to the quadratic level. Again, the number of terms 
is reduced by keeping the sum of the sea quark masses fixed; and we again find the 
number of coefficients reduced by the constraint that the A mass is independent 
of 6fis, and the fl mass independent of Sfii. The most general quadratic formula, 
consistent with the above constraint, and with SU{3) symmetry, is 



Ma = Mo + 3A6ni + Bo6mf + 3Bi6iJ,t 
M5 
M 



= Mo + A{26iJi + djjs) + Bo6mf + 5i(25/if + 5/i^) + 52(5/i, - 5/i/)^ 
= Mo + A{Sfxi + 2(5/i,) + BoSmf + Bi{Sfx^ + 25 ^i]) + B2{5^is - Sfxif 
Mn = Mo + 3ASfXs + BoSmf + 3BiSfil. (53) 



These formulae apply when the sum of the sea quark masses is held constant, 
^{2mi + nis) = mo, but the valence quark masses are completely free, because at 
this level (terms up to second order) a restriction on valence masses would not 
lead to any reduction in the number of free parameters. 

We can check the formulae flS^ by forming the SU (3) mass combinations of 
Table [TTl and checking that in each case only polynomials of the appropriate 
symmetry appear in the answer. 
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There are some combinations of the partially quenched masses, eqs. (153|1 . 
which have simpler dependences on the valence quark masses and only involve a 
small number of fit parameters. Examples include 

- Ma + Ms* + Mh. - Mf^ = 2S2(5/x, - 

Mh. - Ms. = A{5fXs - 6 111) + Bi{5iil - 5fi^) 
Mh-Ma = 3A{6fis-Sfii) + 3B,{Sfil-6fi^). (54) 

Let us use these formulae to illustrate how partially quenched measurements 
might help us fit masses on the constant sea-quark mass line. If we only use 
unitary data, with i(2mf + mf) = m^, we are limited to points on the line 
between the endpoints {mf^mf) = (0, 3m^) and (mf,m^) = (|m^,0), which 
means that the quark mass splitting is limited to the range — < (mf—mf) < 
StUq. In the partially quenched case, we can increase /i^ without having to 
decrease fii, so we can make the splitting fj,f — fif much larger than — mf , 
which gives a much better lever-arm to determine Bi and B2. 

Because the decuplet baryon has a particularly high degree of symmetry we 
can give an alternative derivation of eq. fl53|l . with less explicit reference to the 
flavour group. Consider a decuplet baryon made from quarks of type a, b and c. 
Since the decuplet is a fully symmetric representation, the mass of the baryon 
must be a symmetric function of Sfia, ^fJ'b and 6fic- 

The mass can also depend on the sea quark masses Smu,Sm(i,Sms, but this 
dependence must be in a flavour singlet way. If we keep the sum of the sea 
quark masses fixed, the first singlet polynomial allowed (see Table [3]) is liSm"^ + 
6m'^ + (5mf), where we have chosen the prefactor | to lead to a tidy expression 
in the limit when = rrid- Since this mass polynomial is flavour singlet, it can 
only appear in the coefficient of the identity matrix - i.e. it must make the same 
contribution to every decuplet baryon. We therefore know that the mass formula 
for the abc decuplet baryon must have the form 

M{ahc) = Mo + BoKSml + dm^ + 5ml) + FsymiS^a, S^t, SjJc) 

Mo + BoSmf + FsymiSiJ^a, Sub, ^/^c) , (55) 

in the 2 + 1 limit. 

Now let us consider what terms are possible in the completely symmetric 
function Fgym- There is only one symmetric linear term, 6fia + SfJ'b + ^f^c- There 
are two symmetric quadratic terms we can write down, Sfi"^ + Sfif + and 
S^aSfi'b + SfiaSfic + SfibSfic, OY linear combinations of these two terms. We get 
simpler final expressions if we choose the basis 6fil + 6fil + Sf^l and Sfil + 6fil + 
Sfj^l — SfiaSfib — Sfia^fJ'c — SfibSfic, giving us the general mass formula 

M{abc) = Mo + A{Sfia + 5fXb + 5fXc) + Bo5mf + Bi {Sfil + 5^1 + 5/i^) 

+B2{5^il + 5^il + 5^il - dj^adfib - Sjj^aSfic - 5^ib5^c) ■ (56) 
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Calculating M{lll), M{sll), M{ssl) and M{sss) gives the partially quenched mass 
formulae in eq. (!53|) . Remember that these formulae are only complete for the 
case ruu + rud + rus held constant. If the average sea quark mass is allowed to 
vary, more terms become possible, including 'mixed' polynomials which contain 
both 5mq and d^iq. With fixed average sea quark mass, such mixed polynomials 
do not arise until the cubic order. 

This argument gives us the same result as the full group-theory derivation, 
though in some sense it explains less. For instance, it does not immediately 
explain why the mass combinations of Table [TT] give particularly tidy polynomials, 
or where the hierarchy in eq. ( H2l) comes from. 

4.2 PQ octet mesons 

As in the previous section, we find the meson mass formula by constructing the 
most general matrix consistent with 5'f/(3) symmetry, and the constraint that 
the partially quenched pion mass must know nothing about yU^. The resulting 
mass formulae for partially quenched mesons take the form: 

Ml + 2a5/i/ + l3o5mj + 2/3i5/if 

Ml + 2a5fXs + Mm] + 2(3i5yil . (57) 

Again the rjs is the meson made of a partially quenched s^a/s^a/ quarks (i.e. the 
'strange pion') which in the partially quenched framework can be observed and 
can yield useful information about the extrapolation constants. The PQ M^^ can 
thus be obtained from the PQ by simply changing /i; — )■ /x^ which changes 
the top row of eq. (!57|) into the bottom row. 

Some useful combinations, which avoid the delicate r] sector, are 

Ml -Ml = a{6fis-6fii) + PiiSfil- 6 P2{6^^s-Si^i)'' 
2Ml + Ml = 3Ml + a{4:5fXi + 25fXs) + ^(3o5mf 

+/3i(45/i2 + 25fi1) + 2^2(6 fis - 5/^0' . (58) 

Ml — Ml is useful as a measure of the quark mass splitting, 2Ml + Ml as a 
quantity which is nearly constant along our trajectory. 

The same form, mutatis mutandis, applies to the other meson octets, e.g. the 
p, K*, (j) system. We thus have 

Mp = Mo + 2a6fii + f3oSm^ + 2/3i5/if 
Mk* = Mo + a{6pi + 6ps) + (3oSm^ + (3i{SiJ,f + 6pl)+f32{Sps-Spi)^ 
M^^ = Mo + 2a6ps + Mmf + 2/3i6fil, (59) 

following the pattern of eq. (1571) . 



M: 



K 



Mi 
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We can give a similar elementary argument to derive the partially quenched 
mass formula for mesons with different quarks (e.g. the K* and the charged vr 
and p). Consider an ah meson. It must have the same mass as its antiparticle, the 
ha meson. So, by the same argument as in eq. (l55l) we will have a mass formula 

M{ah) = M{ha) = Mq + PoliSml + 6ml + ^^1) + Fsym{Sfia, Sfib) . (60) 

The only linear symmetric polynomial is + 6fif„ while there are two indepen- 
dent quadratic possibilities, which can be chosen to be Sfi^ + Sul and {6fj,a — Sfj,b)'^. 
With this choice we get the meson mass formula 

M{ab) = Mo+a{6iJa+SiJb)+/3ol{Sml+6ml+6ml)+(3i{6fil+6ijl)+/32{6iJa-Sfiby , 

(61) 

which reduces to eqs. f l57|) and fl59|) . (For pseudoscalar mesons we expect to 
have a smoother Taylor expansion than M itself, so we keep the form of eq. fl^T]) . 
but apply it to the squares of the pseudoscalar masses.) This formula will also 
apply to the ss meson with annihilation 'switched off', i.e. with disconnected 
diagrams dropped. 

We cannot see a way to extend this argument to include the ?78 or 08 mesons. 
Using our full group argument we find the formulae eq. fH^ and eq. (Hlj) . which 
involve an additional quadratic parameter, fS^. Physically it is very reasonable 
that the rjg and 08 should have a term that cannot be linked by symmetry to 
the mass of a meson with two different valence quarks — the rjs and 08 have 
contributions from qq annihilation, which is absent in the other mesons of the 
multiplet, so it would be surprising if symmetry could completely determine the 
masses of these 'central' mesons. 

4.3 PQ octet baryons 

The number of free coefficients in the meson case was reduced by the requirement 
that K and K have the same masses. However, there is no similar constraint 
linking N and H, so more coefficients are allowed, both at the linear and quadratic 
levels. Arguing as before, we find 

Mn = Mo + 3Ai6fii + Bo6m^ + 3Bi6fif 

Ma = Mo + Ai{26fii + 6i2s)-A2{6i2s-6i2i) + Bo6m^ 

+B^{26fi^ + 6fil) - B2{6i4 - 6fi^) + B^idfi, - 5^iif 
Ms = Mo + Ai(25/ij + 5/i,) + A2(5/i, -5/i/) + 5o5mf 

+5i(25/i2 + 5iil) + B2{5iil - 5^f) + 53(5/^, - 5iiif 
Me = Mo + Ai{5fii + 25i2s)-A2i6fis-6fii) + Bo5mf 

+B,{6fi^ + 26fil) - B2{6fil - 5/1^) + Bs{6fi, - 6fiif . (62) 
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As usual, the nucleon mass has been made independent of 5^s- Some useful 
combinations, which only depend on a few parameters, are 

2M^ - Ms - 3Ma + 2Mh = (^3 - 354)((5;U, - 5^/)^ (63) 
Ms-Mj, = (Ai-2A2)((5/i,-5/i/) + (5i-252)(5/i^-5/i?). 

As mentioned previously, we can check that when fi ^ m (i.e. return to the 
'unitary line') then these results return to the previous results of eqs. (143!) - fH6|) . 

As with the decuplet baryons, and the partially quenched mesons, we can 
give an alternative derivation, with less use of explicit group theory. However, 
the argument for partially quenched octet baryons is slightly more complicated, 
because there are fewer symmetry constraints. We will first consider the baryons 
of the type aab, two valence quarks of flavour a, and one of flavour b. These 
form the outer hexagon of the octet diagram. As before, the sea quarks must 
contribute equally to all masses in the octet, the mass formula must take the 
form 

M{aab) = Mq + Bol{6ml + 6ml + 5ml) + F((5/ia, 5/if,) , (64) 

but now the dependence on the valence quark masses is not symmetric under 
a ^ 6, so the function F need not be symmetric. This means there are two 
independent linear terms (we choose {26fia + Sfib) and {6fih — Sfia))- There are 
three independent quadratic terms, we choose (25/i^ + 5/i^) and (5/i^ — 5/i^) (to 
mirror the pattern of the linear terms) and {Siib—SnaY- Thus, the general formula 
for the aab baryons can be written 

M{aab) = Mq + Ai{2SfXa + 5fXb) + A2i5fXb - Sfia) + Bol{Sml + Smj + Sm^) 

+B,{26fil + 6fil) + B2{6fil - 6fil) + B^{5^ib - Sfiaf ■ (65) 

Taking the cases M{Ul), M{lls) and M{ssl) gives the A^, E and H masses of 
eq. ([62]). 

We have not found an equally simple argument for the mass of the A. The 
result of the group theoretical calculation, as set out in Appendix |Bl is 

Ma = Mo + Ai{26 1^1 + 6 12s)-A2{6i2s- 6121) + Bol{6ml + 6ml + 6ml) 

+Bi{26fi^ + 6 111) - B2{6iil - 6fi]) + B^{6fis - 6fii)^ . (66) 

Most terms are related to terms in the S mass, the Mq, Ai, Bq and Bi terms are 
the same for A and S, while the A2 and B2 terms have opposite signs for the A 
and S. However, for the term {6fis — 6^1)'^ there is no connection between the 
coefficient in the A mass and the coefficient of this term in the other masses, and 
we need to introduce a new parameter, that can only be determined by simulating 
the A. We can understand this from Table [HI There is a particular combination 
of singlet, 8^ and 27-plet matrices which gives for all the baryons in the outer 
ring, and only acts on the central baryons. Clearly, the coefficient of this matrix 
only appears in the A mass, so at the quadratic level, we can no longer predict the 
A mass from the other masses in the octet. We have a similar situation with the 
mesons - there is a quadratic coefficient that we only see in the r^g mass formula. 
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4.4 Generalising a constant m formula 

We have stressed the advantages of keeping the average sea-quark mass, m, con- 
stant when approaching the physical point. This leads to simpler extrapolation 
formulae, and results closer to the physical values for flavour singlets and for par- 
tially quenched calculations. If we want to move away from the surface m = mo, 
for example to consider a completely different trajectory, such as = constant 
or mf = constant, then it would be useful to know how to generalise our constant 
m formulae to cover the full parameter space. 

The procedure is simple; every constant parameter in our formulae becomes 
a function of m, which we can then Taylor expand around m = mo. Taking as 
a flrst example our cubic formula eq. (fT6|) for a flavour singlet quantity (such as 

— = a + /3{Sm1 + Sm^ + Sml) + 7 SmuSmdSrUs , 
a 

a + aim — mo) H — -a"{m — mo)^ H — -a"'{m — mo)^ 
2 . 3 • 

+(3{5m\ + + 5ml) + (3'(m — mQ){5ml^ + 5m\ + 5m1) 

+7 5mu5md5ma , (67) 

yielding a cubic formula with 7 parameters. The extra polynomials appearing in 
this formula are the 'unticked' Ai polynomials in Table [31 

We can take a slightly more complicated example, the partially quenched 
formula for the S baryon. In eq. ( l62l) we give the quadratic formula, valid with 
m held constant. Thus to construct the quadratic formula without this constraint, 
we must Taylor expand the parameters Mq, Ai,A2, giving 

Ms = Mo + M^(m-mo) + ^M^'(m-mo)^ + Ai(2(5/i, + (5/i,) 

+A2{Sfj,s - Sfii) + A'i(m- mo) (2(5/1/ + + A'^im - mo){6fis - Sfii) 
+BoSm^ + 5i(25/i2 + 6fil) + B2{5nl - 6fi^) + 3^(6 fis - Sfnf . (68) 

Note that this formula contains 'mixed' polynomials such as {m — mo){6fis — Sfii), 
involving both valence and sea quarks. If we restrict ourselves to the constant m 
surface such mixed polynomials only show up at the cubic level. 

4.5 The usefulness of PQ 

There are several possible advantages to considering PQ results. 

1. The coefficients that appear in the expansions about the flavour symmetric 
line in the PQ case are the same as those that appear on the 'unitary' case. 
Hence this may be a computationally cheaper way of obtaining them. 
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2. PQ results can be helpful in choosing the next point to simulate, because 
the meson masses measured in the partially quenched approximation are 
very close to those found in a full calculation, giving us a preview of results 
on the next simulation point. We understand theoretically why this works 
well on our trajectory, with fn held fixed. The reason is that the effect 
on the sea of making the u and d quarks lighter is largely cancelled by 
the effect of making the s quark heavier (the cancellation is perfect at the 
flavour symmetric point). Therefore partial quenching works best when 
only the non-singlet part of the quark mass matrix is varied (as is the case 
here). If the singlet part (the average sea-quark mass) is changed, there is 
no compensation, and the partially quenched results are less reliable. 

3. We can use partial quenching to get a good estimate of results at the phys- 
ical point, by taking configurations generated with quark masses some dis- 
tance short of the physical point, and then at the measurement stage using 
valence quarks chosen to give the physical vr and K masses. Important 
physical effects, such as the light pion cloud, would be incorporated in the 
results. The effects of partially quenching can be further reduced by repeat- 
ing the calculation with several choices of sea-quark masses, and making an 
extrapolation towards the physical sea-quark mass values. 

4. It is necessary in the determination of strange (or ss) mesons without dis- 
connected pieces. 

We can also show that on this trajectory the errors of the partially quenched 
approximation are much smaller than on other trajectories. In leading order xPT 
(terms linear in the quark mass), the suggested procedure (valence quarks at the 
physical value, sea quarks anywhere on the physical constant m^ea) is exact. See 
Table VIII of |2lj for the leading order formulae for both octet and decuplet 
baryon masses. At this order partial quenching moves all the octet baryons by 
the same amount, and all decuplets also move together. The leading order partial 
quenching errors are 

MZ-M^ec = -6aM(m-m*), (69) 

(using the notation of [21] for quark masses and the a coefficients). The super- 
script * denotes quantities at the physical point. Since we have tuned m to be 
equal to the physical value, the partial quenching error vanishes on our trajectory, 
but not on other trajectories, which vary m. 

We can give a partial derivative argument, like that of section [1] or [Ij, which 
explains why this is so. Take the proton mass as an example, but any quantity 
will work the same way. The proton mass will depend on the valence quark 
masses and the sea quark masses, so we can write 

M^{^u,fJ'd;mu,md,ms) . (70) 
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is the mass of a partially quenched hadron calculated on a sea background. 
The dependence on the three sea masses must be completely symmetrical, unlike 
the dependence on valence masses. At the symmetric point 

dM^ dM^ dM^ 
druu drud dnxg 

so if the sea quark masses are changed in a way which preserves m, while the 
valence masses are held constant, will not change (to leading order). We 
can see these benefits of the constant m procedure in the mass formulae of this 
section. If we make our valence quark masses equal to the quark masses at the 
physical point, the only difference between the partially quenched hadron mass 
and the physical hadron mass comes from the or /5o term in eqs. (l53l) . (1571) . 
(1621) which gives a quadratic mass shift 

M*-M^ = 5o[(5mn'-5mH, (72) 

(where Sm] refers to 5mi at the physical sea point), which is one power higher 
in the quark mass than the usual result on other trajectories, eq. fl69l) . This 
partial quenching shift is the same for all particles in a multiplet, so splittings 
are unaffected by partial quenching at this order - we would have to expand to 
cubic terms to see partial quenching errors in the splittings. 



5 Applications to chiral perturbation theory 

Almost all LO (i.e. leading order or zero loop) chiral perturbation theory (xPT) 
results follow simply from flavour blindness, without any input from chiral sym- 
metry. The linear terms in m^, which are usually called LO xPT results, were 
originally discovered by Gell-Mann and Okubo [3 [8], using the (non-chiral) SU (3) 
argument we are using in this article. 

The only case where we need to invoke chiral symmetry is for the pseudoscalar 
meson mass formula, where it is chiral symmetry which tells us that we have 
massless Goldstone bosons if 2 or more quark masses vanish. 

Beyond leading order we cannot derive the xPT result in full solely from 
flavour blindness, but we can still make useful statements about the form that 
higher order contributions must take. 

3/2 

5.1 Decuplet baryon masses at 0(mq ) 

0{my^) xPT is based on one- loop graphs, all with a pseudo- Goldstone boson. So 
we should expect that the individual terms in the xPT answer will be functions 
of or of Mk or of M^g, with no mixed terms (such as M^M|^), which can 
only arise at the two-loop level. 
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As an example, let us examine the 2 + 1 next to leading order (NLO) results 
for the decuplet baryon masses, [22j. Taking the formulae for individual masses, 
and grouping them into the multiplets of Table [HI we know that in each case 
we are only allowed chiral perturbation theory expressions in the corresponding 
multiplet: 

4Ma + 3M2. + 2Mh* + 



-2Ma + Mh* + Mq 



4Ma - 5Ms. - 2Mh. + 3Mj 



-Ma + 3Ms* - 3Mh* + 

using the notation of [22] (In particular J^-{Mi) is short-hand for the function 
J^{Mi, — A,/i) defined there). The coefficients on the right-hand side of eq. (|76ll 
must follow the pattern of Table [121 but any function of the meson masses is 
allowed. We proved a weaker version of this result in [1], using the permutation 
group instead of full SU{?>). 

The meson mass matrix, 8 (g) 8, contains no 64, there are no possible 1-loop 
terms to place on the right-hand side of eq. f[7^ . so this mass combination must 
be zero in NLO xPT. We have already noted in Table [11] that this combination 
has a Taylor expansion beginning at 0{5m^) and is very small experimentally, 
eq. (|43D. 
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(75) 



-J-_(M^) + 4J-_(Mk) - 3J-_(A/r,J] 

(76) 



5.2 Relationships between expansion coefficients 



We now investigate the relation between the parameters of xPT and the Taylor 
coefficients in our approach, eqs. (H3|) - ( H6|) . 

For example for the pseudoscalar octet, using the 2 + 1 results in [23] and 
assuming their validity up to the point on the flavour symmetric line (the kaon 
mass is always smaller here than the physical kaon mass), we flnd 



Ml 



X 



16x 
Jo 



(3L4 + L5- 6L6 - 2L8) + 



24^ "^Af 



38 



a 



= Qo 



^ (3L4 + 2L5 - 6L6 - 4L8) 



^0 = 




67r2/o' 




Po 


-32 (L5 




Ql 

f2 


16 (L5 - 



(77) 



where Xq = '^Qof^q with Qq = B^Z^ so that here we have x = 2Qo(l + ctz)^ 
which is kept constant. The LjS are appropriate low energy constants or LECs. 

We first note that when expanding the xPT about a point on the SU{3) 
flavour symmetry hue gives to leading order only one parameter, a as expected. 
(This means, in particular, that the flavour singlet combination, X^, vanishes to 
leading order, as discussed previously.) Secondly, while we can fit to a and f3Q, 
Pi and /32, it will be difficult to determine the individual LECs. The best we can 
probably hope for are these combinations. 



5.3 Chiral non-analytic behaviour 

We brieffy discuss the question of how chiral logs, or other chiral singularities, fit 
with the Taylor expansion. The answer is that the chiral singularity should show 
up in the large-n behaviour of the coefficients of 5m^. 

For example if we make a Taylor expansion about x of the singular term 
X/^ln(x//A^) which occurs in the chiral expansion of M^, we would find 

In X. = [X + 5xif I^X + X5XI + bxl+t E , ^^"^ C-^) ^ , (78) 

where Sxi — Xi ~X- So at large n the coefficients of ^xF proportional to 
X^~"/n^. If we look at the first singular term in the baryon mass formula, x/"^ ^ 
and expand, we would get a series ~ x^/'^~^5x^ /n^^"^ at large n. 

This is general, the power of n with which the terms drop off tells us the 
chiral singularity. If the singularity is xnn(xi/A^), with integer p, then the 
series drops hke l/vP'^^. If the singularity is Xp with non- integer q, then the 
series drops like l/n'^+^. (If we have a singularity in ln(x,7/A^) or ln(xAr/A^), 
where Xv = ixi + 2Xs)/3, Xk = ixi + Xs)/2 we just change xi to Xv o^' Xk in the 
Taylor series.) Needless to say, this large-n behaviour would prove difficult to see 
in practice, because the coefficients are small these terms only become important 
close to the chiral limit. 
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6 The path to the physical point 



In section [T] the proposed path to the physical point was introduced. We shall 
now discuss this a little further. 

For the simulation it is easiest to keep the (bare) singlet quark mass fixed, 



m = i(2m; + rris) = niQ = constant 



(79) 



starting from some reference point (m;, m^) = (mo, mo) on the fiavour symmetric 
line. We can use the singlet combinations from Table [10] to locate the starting 
point of our path to the physical point by fixing a dimensionless ratio such as 



4 

Xn 



physical value 



10) 



Note also that at the fiavour symmetric point = M.^\q so this determines our 
starting pion mass (from the experimental values given later in Table [T6l) to be 
~ 410 MeV. 

However the equivalence of eqs. (17^ . (ISUi) is only strictly true at lowest order. 
While at this order it does not matter whether we kept the quark mass singlet 
constant, eq. f l79|) . or a particle mass singlet constant, eq. flHOj) . higher order 
terms mean that it now does. If we make different choices of the quantity we 
keep constant at the experimentally measured physical value, for example 
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Xn 



XI 



XI 
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we get slightly different trajectories. The different trajectories begin at slightly 
different points along the ffavour S'f/(3) symmetric line. Initially they are all 
parallel with slope —2, but away from the symmetry line they can curve, but will 
all meet at the physical point. (Numerically we shall later see that this seems to 
be a small effect.) 

An additional effect comes from the choice of Wilson lattice fermions. The 
physical domain is defined by 



which using eq. ([5]) 

mi > 



translates to 
1 . 



m, 



m„ 



m.. 



> 
> 








m., > 



-nil 



12) 



(83) 



leading to a non-rectangular region and possibly negative bare quark mass. (These 
features disappear of course for chiral fermions when az = 0.) 
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Figure 9: Sketch of some possible paths (red lines) in the mi~ms plane to the physical 
point (m^ , m* ) . Due to eq. ([83]) the mf and axes are not orthogonal when plotted 
in the bare quark mass plane. 

These two features are sketched in Fig. [9l which shows possible paths in 
the mi-ms plane starting from the flavour symmetric line. In this figure the 
horizontal and vertical axes are the bare quark masses, mi and respectively. 
Because of renormalisation effects, eq. fl83|) . lines of constant renormalised mass, 
mf = constant or mf = constant, will be at an angle. This is in contrast with 
Fig. [H where renormalised masses were used as the axes. 



After the general discussion of SU (3) flavour expansions described in sections |2H5] 
(which are lattice independent), we now turn to more specific lattice considera- 
tions. 

7.1 Lattice simulations 

We use a clover action for 2 + 1 flavours with a single step of mild stout smearing 
as described in Appendix Further details are given in [21] together with 
a non-perturbative determination of the improvement coefficient for the clover 
term, using the Schrodinger functional method. 
The bare quark masses are defined as 



7 The lattice — generalities 
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where vanishing of the quark mass along the SU{3) flavour symmetric line deter- 
mines kq-c- We then keep m = constant = mo which gives 

So once we decide on a this then determines Kg. 

How accurately must we satisfy eq. (185!) ? In choosing suitable {ki,Ks), the 
natural scale is to say that changes in m should be small when compared to fn 
itself, i.e. \fn — mo\ -C mo which gives 



6 \Ki Ks Ko 



«--— , (86) 



which is satisfied if we give all our ks to 6 significant figures. 

Furthermore note that we are not expanding about the chiral limit, but have 
expansions around a flavour symmetric point which does not require knowledge 
of kqjc- This follows as 



7) 



bvTiq = ITLq — m 

- -{--- 

2a \Kq Kq 

HMC (hybrid Monte Carlo) and RHMC (rational HMC) were used for the 
2 and 1 fermion flavours respectively, [25], to generate the gauge configurations. 
We note the following in connection with the simulations and our path choice: 

• The simulations should equilibrate quickly from one point to another along 
this path, because the effects of making the strange quark mass heavier 
tend to be cancelled by making the u and d quarks lighter. 

• The simulation cost change should be moderate for this path. This may be 
motivated by the following crude cost argument. Modelling the cost, C, as 

Coc^ + -^, (88) 



where k is the relative cost of the two algorithms, gives on the line afn = 
constant 

^(^) ^ n M + \ 0, ' (89) 

{l + az)+i (l + az)-2^ 

with 

C = Smi/m, (90) 

(alternatively we could consider M^/X^). The cost C{^)/C{0) from the 
symmetric point ,^ = 0, is plotted in Fig. [TUJ There is little change in a 
reasonably large range of ^ starting from ^ = 0. 

Both these points are indeed found in practice (at least approximately). 
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Figure 10: Simulation cost C(^)/C(0) against ^ where ^ = 6mi/rn (with, for example 
az = 0.75, cf. eq. (|104p . and k = 1). The symmetric point is denoted by a (red) filled 
circle. Very roughly we need to reach a region where ^ < —1.5. 

7.2 0(a) improvement of the coupling constant 

0{a) improvement leads to a change in the coupling constant via [26] , 



Wilson glue and a = (i.e. no-stout smeared links) Wilson-Dirac fermions the 
lowest order perturbative result is hg = 0.01200n/5fQ -|- O^Qq), [26], which is small 
but increasing with Uf (here n/ = 3). A crude estimate was made in [27] and 
indicated a possible 1-2% effect (but with considerable uncertainty). 

In general from eq. (1911) when we vary a quark mass then Qq must be changed 
to keep go constant. However for our choice of path (m = constant) the relation 
between go and go is fixed, so only a small overall shift of results might be neces- 
sary - nothing else changes as we traverse our path. Therefore in the following 
we shall not consider the effect of bg any further. 

7.3 Hadron 'sources' and 'sinks' 

The operators (or interpolators) used to determine the hadron masses are uni- 
formly taken to be Jacobi smeared ([28] and [29] (Appendix C)) and to be non- 
relativistic, NR, ([30], [31] and [32] (Appendix C)). Specifically, we consider the 
following hadron sources and sinks: 

• Pseudoscalar meson octet 




(i(x)75'u(x) 
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Mk{x) = s{x)-f5u{x) 



(92) 



• Vector meson octet 



Mpi{x) = d{x)-fiu{x) 
MK*i{x) = s{x)jiu{x) 

M^^i{x) = s{x)-fis{x) (93) 



Baryon octet 



BNa{x) = e'^'^ul{x)[u\xf''C^,d%x)] 
Ba^{x) = e<^^'^(2<(a;) [M''(a;)^-C75d^(x)] 

+C(a:) [u\xf''C'y,s''{x)]-uUx) [d\xf^C^5s'{x)]) 
i3sa(x) = e^'-Xl^) h'W^"C'75S^(:r)] 

Bsa{x) = e"''^!^) [5^x)^-C75«^(a;)] (94) 

• Baryon decuplet 

BAaix) = e"'^(2<(x) [M^(x)^^C7_rf^(x)] +rf:^(a;) [M'(a;)^^C7_M^(a;)]) 

Bj:*aix) = e"'^(2<(a;) [M^(a;)^^C7_s^(x)] + <(x) [M''(a;)^°C7_M'=(a;)] ) 

BE^aix) = e"'"(2s:i(x) [s^(a;)^°C7-ix^(a;)] +<(a:) [s\xf''C'y.s''{x)]) 

Baa{x) = e'''''sl{x)[s''{xf''Cj_s%x)] (95) 

where C = 7274 and 7_ = §(72 + ^7i) and the superscript is a transpose in 
Dirac space. The u and d quarks are considered as distinct, but of degenerate 
mass. 

The correlation functions (on a lattice of temporal extension T) we use are 
given by 

\ y X I 



3^3 

y 



\ ? S I 
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\ y X I 

^ \ y X I 

oc Ae-^^o* + Be'^^'^o^^''^ , Aq = A, E*, S*, ^] (96) 

with Tunpoi = |(1 + 74) and Tpoi = T^npoii^ + HsTs)- M' is the lowest excited 
state with opposite parity to M. 



8 The lattice — results 

All the results given in this article will be at /3 = IO/qq = 5.50, a = 0.1, together 
with Csw = 2.65, see Appendix [Dl (This (3 value was located by an initial series of 
short degenerate quark mass runs, to give a rough idea of the associated scale.) 
The hadron masses will be given below in a series of Tables [T9] - [30l The runs on 
24^ X 48 lattices have 0(2000) trajectories, while the runs on 32^ x 64 lattices 
have O(1500)-O(2000) trajectories. Errors are determined using the bootstrap 
method. Experimental values of the hadron masses are given in section O 



8.1 Locating kq and the nig— mf plane 

From the discussion in section [6] for our path choice, we must first determine the 
starting value on the flavour symmetric line. A series of runs along the SU{3) 
flavour line determines this point, kq, by checking when X^/Xg, S = N, A, p 
are equal to their physical values, see eqs. flHU]) . fIST]) . (This would also include 
5" = r if we have previously determined the physical value of tq.) On the flavour 
symmetric line obviously all the particles in the multiplet are mass degenerate, 
so for example taking S = N means that 



{aM^f XI ' ^^^^ 

where, to emphasise that the left-hand side are the lattice measurements, we 
temporarily include the lattice spacing. (Again, the star denotes the experimental 
value.) 

Once we have located a promising (or better a small range of kqs) then we 
keep m = constant and pick appropriate {ki,Ks) values, using eq. ( I85l) . Again, 
setting X^/Xlf = physical value, eq. ( IHOl) can be re- written as 



2Ml-Ml_ Ml , _.Xl 

- Cat - l-TT ' Cn — 6 — 2- 



(9^ 
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considering for the present only lowest order in the flavour expansion. In Fig. [TT] 
we plot (2M|^ — M^)/Xjf versus M^/Xjf. This is equivalent to plotting mf 
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Figure 11: {2M'^-M^)/Xjf versus M^/Xj^ for kq = 0.12090 (left panel) and 0.12092, 
0.12095 (right panel). The dashed black line, y = x represents the SU{3) flavour 
symmetric line. Filled points are on 32^ X 64 lattices while open points are on a 24^ X 48 
sized lattice. Shown are points on the flavour symmetric line (coloured orange) followed 
by results with rn = constant (coloured violet, left panel; blue and green, right panel). 
The fits are from eq. ()98p with cat a free parameter. The physical value is denoted by 
a (red) star. 

against mf because from LO xPT, oc mf and 2M|^ — oc mf. 

Note that simulations with a 'light' strange quark mass and heavy 'light' 
quark mass are possible - here the right most points in Fig. [11] In this inverted 
strange world we would expect the weak interaction decays p — > S or A. 

Also shown in Fig. [11] are fits using eq. ( p8|l leaving free parameter 

starting from the flavour symmetric points 

Ko = 0.12090, Ko = 0.12092, fi:o = 0.12095, (99) 

(the latter two points are reference points). It is seen from the figure that this 
range covers the possible paths to the physical point. There are two observations 
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to be made. Firstly we note that there does not seem to be much non-hnearity in 
the data, i.e. the leading order in the expansion about the flavour symmetric line 
already seems sufficient. So if c^v = 3(X^/X^)|* then the lines would go exactly 
through the physical point. Also this means from the discussion in section [6] that 
using other singlet scale choices should lead to a similar result. Secondly, as noted 
before at the end of section |3l as the expansions have been derived using only 
group theoretic arguments, they will be valid for results derived on any lattice 
volume (though the coefficients of the expansion are still functions of the volume). 
So here, to test this, we have made separate fits for the two volumes — 24^ x 48 
and 32^ x 64. Indeed this shows that finite size effects are present but small. 

Thus our present conclusion is that k in the range 0.12090 - 0.12095 is within 
a few percent of the reference kq. Most of the results reported here will be at 
Ko = 0.12090. 



8.2 Determination of kq-c, az 

Although not strictly necessary, we briefiy indicate here the determination of Kq-c 
and az to illustrate some of the discussion in section |H1 Using lowest order xPT 
(i.e. the fact that the pion mass vanishes if the masses of the light quarks vanish) 
and 

{aMj,Y oc amf oc am; + azofn , (100) 



JVS 
m 5 



where the constant of proportionality from eq. (1771) is 2aa = 2aQo = ^aB^Z, 
we first determine Kq-c (the critical hopping parameter on the fiavour symmetric 
line) as defined in eq. (^^. In Fig. [12] we show the plot of (aM^)^ versus 
together with linear fits. For the fiavour symmetric case, from the blue points we 
find 

— = 8.25768(23) , or kq.c = 0.121099(4) , (101) 

which is in good agreement with the Schrodinger functional determination, see 
Appendix [Dl Note that for ki < kq-c the bare arriq is negative, but the renor- 
malised is always positive, eq. fl83|) and Fig. [91 This occurs for the last point 
on the 32^ x 64 line in Fig. [HI 

az can then be estimated using the m = constant line as here (aM,r)^ vanishes 
at Kc giving 

I - ^) 

az = - "'"^_ (102) 

am 11^ 



3 



Using the 32 x 64 data only (green points) gives 

— = 8.24727(17), or = 0.121252(3) . (103) 
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Figure 12: (aM-^)'^ versus 1/ki for both the flavour symmetric case (blue points) and 
keeping m constant (green points). 24^ x 48 volumes are opaque circles and 32^ x 64 
volumes are filled circles. The kq = 0.12090 points are highlighted in red. The chirally 
extrapolated values from the linear fits are denoted by stars. The horizontal (orange) 
filled circle is the Schrodinger functional estimate. 



Hence this gives here 

az~ 0.76. (104) 

Note that the determination is quite sensitive to small changes in kq-^ and Kc- We 
conclude that for clover fermions at our lattice spacing az is indeed non-zero. 

8.3 Singlet quantities and the scale 

We take Fig. [11] as a sign that singlet quantities are very flat and the fluctuations 
are due to low statistics. We now investigate this further. In Fig. [13] we show 
aXs for S = A, N, p and vr against M^/X^ for kq = 0.12090 (left panel) 
and comparison results for 0.12092, 0.12095 (right panel) together with constant 
fits. This indicates that other singlet quantities are also rather flat (we interpret 
variations in Xa to be due to statistical fluctuations). Again fits are made for 
each lattice volume separately. 

We showed in section 12.21 that singlet quantities must have zero derivative at 
the symmetric point. A second derivative would be allowed, but we see from the 
left-hand panel of Fig. [13] that it must be very small. 
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Figure 13: aXs for S = A, N, p and vr versus M^/X^ for kq = 0.12090 (left panel, 
circles) and 0.12092, 0.12095 (right panel, circles and squares respectively) together 
with constant fits. Filled points and lines are for 32'^ x 64 lattices, while opaque points 
and dashed lines are for 24^ x 48 lattices. (In the right panel the lower filled points 
and lines are for kq = 0.12095.) The physical point corresponds to the dashed line 
at Ml/Xl = Ml/Xl\*, while the symmetric point corresponds to M^/X^ = 1 in the 
figure. 



8.3.1 Finite size effects 



In Fig. [13] there are again indications of relatively small finite size effects. We 
now briefly investigate this a little more. While we do not attempt to formally 
derive a formula here, we do have the obvious constraint that the finite size 
Xs must also be fiat at the symmetry point (symmetry arguments apply in any 
volume). The various possibilities are given in eq. pij) . The first /-term counts 
the contributions of the kaons and charged pions; the second (7-term is irrelevant 
because the strange pion is fictitious. The third /i-term accounts for the 1] and 
TT*^. So it is likely that the first term is dominant because there are more particles 
exchanged (the functional forms are all likely to be similar). So if we only want 
a rough estimate (for the x-axis of the plot) then we shall just choose the first 
term. 
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Thus from eq. (HT]) and as we shall consider only the lowest order term from 
eq. ( HHj) , we expect the finite size functional form to be 

XsiL) =Xs{l + cslihiM^) + 2/l(Mx)]) . (105) 

Lowest order xPT, [331 El] indicates that a suitable form for fiiM) is 

/^(M) = (aM)^ 3-, meson, 

(ML) 2 

-ML 

h{M) = {aMfj—— baryon. (106) 
In Fig. [H we plot Ul{.M^) + 2/i(Mi^))/3 against aX^ for 5 = A, A^, p and vr on 
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Figure 14: 0X5 versus (/l(M^) + 2/l(M;^))/3 for = 0.12090, with 5 = A (circles), 
(squares), p (diamonds) and vr (upper triangles). The left-most clusters of points are 
from the 32^ x 64 lattices (L = 32a), then follow 24^ x 48 and finally 16^ x 32 lattices. 
The dashed lines are linear fits. 



32^ X 64, 243 X 48 and 16^ x 32 lattices for = 0.12090. The fits are linear. A 
reasonable agreement is seen. (The noisiest signal is for 5 = A.) We see that the 
extrapolated (i.e. L — )■ 00) results are very close to the largest lattice results (i.e. 
32^ X 64), so we conclude that using the largest lattice size available should only 
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introduce small errors. We shall also go a little further and assume that finite 
size effects for masses are similar to those of Xs for each mass of the appropriate 
multiplet. Thus we shall later consider ratios M/Xs for all the available lattice 
data; finite size results then tend to cancel in the ratio. 

8.3.2 Scale estimation 

One advantage of our method is that Xs remains constant and can be used to 
determine the scale. We do not have to first extrapolate to the physical limit in 
distinction to other methods. 

The result of section 18.3.11 is that the largest volumes available seem to have 
small finite size effects, so we now simply take the largest volume available. In 
Fig.[T5]we plot aXs/Xg, for S = N, A, p, n using the largest volume fitted results 
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Figure 15: aXs/Xg against I/kq for 5 = A (green squares), (red circles), p (violet 
diamonds), vr (blue upper triangles) with kq = 0.12080, 0.12090, 0.12092, 0.12095 and 
0.12099. 

from Fig. [T3] (together with smaller data sets for kq = 0.12080 and kq = 0.12099). 
The experimental values of Xg are given in section [9711 This ratio gives estimates 
for the lattice spacing a for the various scale^. We would expect most variation 

5For example for kq = 0.12090 we find that a = 0.0834(1) fm, 0.0812(6) fm, 0.0795(3) fm 
and 0.0771(3) fm using X^^, Xa, Xn, Xp to set the scale respectively. 
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of the ratio with and convergence to a common scale where the hnes cross, 
assuming 

• the simulation statistics are sufficient 

• all O(a^) corrections are negligible 

• there is little (or no) curvature present in Xs- 

This appears to be the case, with the possible exception of the decuplet scale. 
However this is the channel with the worst signal, and may be showing some 
curvature (we cannot at present say whether there might be large O(a^) effects), 
so presently we just consider the approximate crossing of the other lines giving 
a ~ 0.075 - 0.078 fm. 

As discussed in section 18.3.11 we expect a (partial) cancellation of finite size 
effects (and also statistical fluctuations) within the same multiplet so we shall 
adopt the philosophy when considering the hadron spectrum of ffist finding the 
ratio of the mass to the singlet quantity from the same multiplet. For example, we 
can take as our base singlet quantity as the baryon octet Xn (not only are these 
stable particles under QCD interactions and so might physically be considered a 
good choice, but X^r also has smaller numerical errors than Xj\ on the lattice). 
To translate from one scale to another we then need the ratio aXs/aXjs[. In 
Fig. [16] we plot Xs/Xn for various Xs (with S = A, p, n). Also shown are 
constant fits to the two volumes — 24^ x 48 and 32'^ x 64. The change in the 
ratios between the two volumes is seen to be small. Note also that all ratios are 
close to their physical values. We use the results of the largest volume, which are 
given in Table [151 In the last column of this table we have used the experimental 



Table 15: Lattice ratios of singlet quantities aXs/aX^, S = tt, p, A from 32^ x 64 
lattices. In the last column we have multiplied by the experimental inverse ratio, taken 
from Table [T6l If we had perfect scaling then this ratio should be 1.0. 

values of Xs (as given in Table [TB]) to form the ratio aXs/aX^ x {X]\f /Xs)\* ■ This 
should be one. As can be seen from Fig.[Tni this is the case and Table [T^ conffims 
that the ratios are 1 within a few percent. All this shows that kq = 0.12090 has 
fn very close to the correct physical value. 



Ratio 



Ko = 0.12090 x{Xn/Xs)\* 



CiXj^ j CLX]\f 

aXp/ (iXj\f 
(xX^I (iX]\f 



0.3751(13) 1.049(4) 
0.7200(38) 0.971(5) 
1.219(9) 1.017(8) 
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Figure 16: aXs/aXN versus M^/X^ for S = A (squares), p (diamonds) and tt 
(upper triangles), top to bottom, for kq = 0.12090. The dashed vertical line represents 
the physical value, while the dotted line gives the SU{3) flavour symmetric point. 
Filled points are on 32^ x 64 lattices while open points are on 24'^ x 48 sized lattices. 
Horizontal lines and dashed horizontal lines represent constant fits to either the 32^ x 64 
or 24^ X 48 results respectively. For illustration, we also show the physical values - 
denoted by stars. 



9 Spectrum results for 2 + 1 flavours 

We shall now discuss our lattice results. 



9.1 Experimental values 

As we will compare our lattice results with the experimental results, we first give 
the experimental masses from the Particle Data Group tables [17] . 

To minimise u-d quark mass differences (and also electromagnetic effects) 
for the experimental data, we average the particle masses over isospin I3 (i.e. 
horizontally in Figs. |2l [3]). This gives the experimental values in Table [17] (we 
postpone giving them here in order to display them with the lattice values in 
Table ITTl) . Using these experimental numbers, the experimental values for the 
hadron singlet quantities used here are then given in Table [161 
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Singlet GeV 



XI = 


V(M2 + 2M|,)/3|* 


0.4109 


x; = 


{Mp + 2Mk^)/^\* 


0.8530 


Ajv - 


: {Mn + Ms + Mh)/3|* 


1.1501 


XI = 


(2Ma + Mn)/3|* 


1.3788 



Table 16: Experimental values for the Xs singlet quantities, X*g, S = n, p, N and A. 



9.2 Mass hierarchy 

We now consider the lattice results for the mass spectrum. First we check whether 
there is a strong hierarchy due to the SU (3) flavour symmetry as found in eq. fH2|) . 
namely 

4Ma + 3Ms. + 2Ms* + Mf^ oc (5mi)° singlet 
-2Ma + Ms. +Mn oc 5mi octet 

4Ma - 5Ms* - 2Mh. + 3Mn oc Smf 27 - plet 

-Ma + 3Me* - 3Ms* + oc 5mf 64 - plet . (107) 

In Fig. [T7] we plot these mass combinations (over Xa) against aSmi for hq = 
0.12090. Also shown are the experimental values using the values from Table [T71 
Note the change of scale between the axes. There is reasonable agreement with 
these numbers. Well reproduced, as expected, is the order of magnitude drop in 
the hadron mass contributions with each additional power of Smi. (See [35] for a 
similar investigation of octet baryons.) It is also seen that while (— 2Ma + M=. + 
Mq)/X^ has a linear gradient in drrii, in the other flts any gradient is negligible 
as expected. To check for possible finite size effects we also plot a run at the 
same {ki, Kg) but using a 24^ x 48 lattice rather than 32^ x 64. There is little 
difference and so it appears that considering ratios of quantities within the same 
multiplet leads to (effective) cancellation of finite size effects. 



9.3 'Fan' plots 

We now show a series of plots of the hadron masses from a small quark mass just 
above the flavour symmetric line down to the physical point. As the masses (of a 
particular octet or decuplet) are all degenerate at a point on the flavour symmetric 
line, then we would expect a 'fanning' out of masses from this point. We consider 
second order flts in the quark mass, but show plots using the pseudoscalar mass 
on the X-axis, i.e. from eq. f l43p . Thus we are using the quark mass as an 'internal 
parameter'. As discussed previously at the end of section |2] and in more detail in 
Appendix [C] this is the natural choice. 

In Fig. [18] we begin with the pseudoscalar octet and show / X^ (ttq = vr. 
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Figure 17: (4Ma + 3Ms* + 2Me* + Mn)/XA, (-2Ma + Me* + Mn)/XA, (4Ma - 
5Ms. -2Mh. +3Mn)/XA and (-Ma + SMe- -SMh* +Mq)/Xa (filled circles) against 
6mi together with a constant, linear, quadratic and cubic term in 6mi respectively, as 
given in eq. ()107p . Extrapolated values are shown as opaque circles. Experimental 
values are denoted by stars. The opaque triangle corresponds to a run at the same 
{ki,Ks), but on a 24^ x 48 lattice rather than a 32^ x 64 lattice. The vertical line is 
at the value of 5mi - dm"^ obtained from a quadratic fit to the pseudoscalar octet, 
eqs. ([43l) . (j47|) as described in section [3] and Fig. [TSl 

K, rjs) against M^/Xl together with the combined fit of eqs. fH3|) . fH7|) . A 
typical 'fan' structure is seen with masses radiating from the common point on 
the symmetric line. Note that the right-most point has a small strange quark 
mass and a large 'light' quark mass, so that the order of the meson masses is 
inverted. 

There is however little real content in this plot - the tto = vr line is obviously 
trivial, for the hq = K line the chiral limit and gradient are known as we have 

^^-^ (108) 

(This can also be seen to 0{5mf) by using eq. fj43|) to form Ml/ and M]^/ 
to 0(5m^).) An inspection of Fig. ITS] shows that the numerical results indeed 
follow very well this line, with a gradient of —1/2 and having in the chiral limit 
a value of 3/2. 
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Figure 18: M^^/X^ {tto = vr, K, r]s) against M^/X^ together with the combined fit 
of eq. (j43p for both the 32^ x 64 and 24^ x 48 lattices. The flavour symmetric point 
('sym. pt.') when kq = 0.12090 is denoted as a red point. Experimental values are 
denoted by stars. The opaque triangle corresponds to a run at the same mass but on 
a 24^ X 48 lattice rather than 32^ x 64. 



However the graph does tell us that for the fictitious rj^ particle, there is 
very little curvature which, as this is a constrained fit, must hold for all the 
pseudoscalar octet particles, including the fictitious one. We also note that ratios 
within the same multiplet do indeed tend to give cancellation of finite size effects. 

In Fig. [19] we plot the vector octet multiplet Mp^/Xp against M^/X^ for po = 
p, K*, (pg. Again finite volume effects tend to cancel in the ratio (normalising 
with the singlet quantity from the same octet) and so both volumes have again 
been used in the fit. The combined fit uses eqs. f H4l) . ( H3|) again with the bare 
quark mass being an 'internal' parameter. Some moderate curvature is now seen 
in the extrapolations. Note that as M^^ ^ M^, the physical must indeed be 
almost a perfect ss state, i.e. we almost have 'ideal' mixing. 

Continuing in Fig. [20] we plot the baryon octet M^^/Xjy for Nq = N, A, 
S, S against M^/X^ and similarly in Fig. [2T]we plot the corresponding baryon 
decuplet M^^/X^ for Aq = A, T,* , E* , Q against M^/X^. Although we have 
included quadratic terms in the fit, there is really very little curvature in the 
results. In both these pictures the correct ordering of masses is achieved (also 
the reverse order behind the symmetric point when we have heavy / quark masses 
and light s quark masses). In particular in Fig. [20] we see that the A-S splitting 
is correct. This is a dynamical effect because A and S both have the same quark 
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Figure 19: Mp^/Xp {po = p, K* , (ps) against M^/X^ together with the combined fit 
of eqs. dMl), (03]) (the dashed lines). Same notation as in Fig. 



content. Also in Fig. l2T]Mv* is indeed constant as expected. 

These results show that the Gell-Mann-Okubo relations work all the way from 
the symmetry point to the physical point. 

The masses (using the scale determined by the appropriate Xs) are given 
in Table [171 The results are rather close to their experimental values. (M^^ is 
taken from [36] which uses a quadratic mass formula and ideal mixing which 
is in agreement with the prediction of LO xPT).) However at present we are 
effectively using a different scale, Xs, for each multiplet. If we wish to convert 
these numbers to a base scale, say X^, then they can be converted using 



M. 



aM 



So 



So 



aX 



X X 



N\ 



N 



aM; 



So 



aXs 



X X. 



/ aXs Xn- 
X I X 

\ aXj^ Xs 



(109) 



where the second factor is given in the last column of Table [T^ . These numbers 
are all ~ 1 (within a few percent). However it is to be noted that this causes 
the largest discrepancy to the experimental value. So the largest source of error 



'If using Xn there is an additional factor from the setting of M^/Xl\* on the x-axis 



Ml 
XI 



0.1128 -> 



Ml 

X 



N 



aX 



N 



aXjr 



0.1023. 



The change in the hadron mass due to this is very small (very much smaller than the error 
bar), so we will ignore this here. 
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Figure 20: Mnq/Xn {Nq = N, A, S, E) against M^/X^ together with the combined 
fit of eqs. (l45l) . (l43]) (the dashed lines). Same notation as in Fig. [I8l 
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Figure 21: Maq/^a (Aq = A, S*, E*, Q) against M^/X^ together with the com- 
bined fit of eqs. (H6]) . (^31) . Same notation as in Fig. [181 
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particle 






expt[GeV] 


aM/aXs 


result [GeV] 


= (M^+ + M^o + M^- 


)/3 




0.1380 


0.3359* 


0.145(1) 


Mk = Mk+ = Mk- 




Is 


0.4937 


1.2015* 


0.518(1) 






ss 


~ 0.685 


1.668(3) 


0.720(3) 


Mp = Mp+ = Mp- 




11 


0.7755 


0.9166(73) 


0.759(7) 


Mk* = Mk*+ = Mk*- 




Is 


0.8917 


1.042(4) 


0.863(6) 






ss 


1.0195 


1.184(12) 


0.980(11) 


Mn = {Mp + Mn)/2 




III 


0.9389 


0.8313(77) 


0.956(9) 


Ma 




Us 


1.1157 


0.9621(142) 


1.107(16) 


Ms = (Ms+ + M^o + Ms 


-)/3 


Us 


1.1932 


1.039(5) 


1.195(6) 


Me = (M=o +Mh-)/2 




Iss 


1.3183 


1.130(7) 


1.300(9) 


Ma 




Ul 


1.232 


0.9047(100) 


1.269(17) 


Me* = (Ms*+ + Ms.o + Ms.-)/3 


Us 


1.3846 


1.007(7) 


1.413(14) 


Ms* = (Mh.o + Mh*-)/2 




Iss 


1.5334 


1.102(11) 


1.546(20) 


Mn = M^- 




sss 


1.6725 


1.191(20) 


1.670(31) 



Table 17: The hadron masses. The third column, 'expt', gives the isospin averaged 
masses. (The r/^ mass is taken from [36J.) The fourth column, aM/aXs, gives the 
numerical results from Figs. [T8]-[2TJ (The aM/aX-,^ values for M^^, Mk are exact.) 
The last column, 'result', has used eq. (I109p to convert the scale to the base scale Xat. 

appears to come from the uncertainty in the consistency of different flavour singlet 
quantities used to determine the common scale. 

9.4 Partially quenched results 

We illustrate partial quenching using baryon splittings as an example. The split- 
tings depend mainly on /x^ — /x/ and only weakly (at second order) on other quark 
combinations. In the PQ data shown here, we have points with a large splitting 
between /i^ and reaching up to points where /i^ — /x; is equal to its physical 
value. We can therefore make partially quenched splitting plots reaching down 
to the physical point. 

We have generated partially quenched results on an ensemble with kq = 
0.12090 and lattice volume 24^ x 48. The first baryon octet splitting 'flag' dia- 
gram, Fig. [221 shows just the PQ data. The second. Fig. |23l shows the PQ data 
in grey, compared with the unitary 24^ x 48 data in colour. While only to be 
taken as an illustration, it shows that the PQ data has the potential to be a good 
predictor of real data. 

Of course partially quenched data is not a complete substitute for simulations 
at the physical point, even for splittings. If we take the S-iV splitting as an 
example, we find from eq. fl62|) . 

Mn-Mj, = {A^ + A2){6fii-6iJs) 
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Figure 22: Partially quenched data, (Mat^ - My)/Xn versus (M^ - M|,)/X^. The 
experimental points are denoted by red stars. 
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Figure 23: A comparison between partially quenched and full data from 24^ x 48 
lattices, {Mnq - Mj])/Xn versus (M^ - Mj^)/X'^. Same notation as in Fig. [22l 
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= {Ai + A2){6fii-6fis) 

+{Bi + B2){5^ii - 6fis){5f^i + Sfis) - Bs{5fii - ^/i,)^ . (110) 

If we plot this baryon splitting against the quark mass splitting — 6 fig) the 
Ai, A2 and B3 terms give a simple parabola (as does the B4 term if we consider 
a splitting involving the A). However the Bi and B2 terms do not depend purely 
on the splitting (5/1; — 6^s), they also depend on the quark-mass sum (5/i; + 6fis)- 
Thus the Bi and B2 terms lead to a broadening of data bands in Fig. [22] (two data 
points with the same value of {6^1 — 6^s) may have differing values of {6fii + 6fis)), 
and can lead to the partially quenched data missing the physical point slightly. 
Although we reach splittings {6fii — 6 fig) equal to and even a little larger than 
the physical quark mass splitting, we do this with our light valence quark still 
noticeably heavier than the real u and d quarks, so {6fii + 6fis) > {Sm^ + (5m*) 
at our end-point. The above argument still applies (with minor modifications) if 
we use the partially quenched meson mass difference — Mj. as a substitute 
for {6^1 — 6^s) on the figure's x-axis. 

10 Conclusions 

We have outlined a programme to systematically approach the physical point in 
simulations of QCD with three flavours starting from a point on the SU (3) flavour 
symmetric line by keeping the singlet quark mass constant. As we move from the 
symmetric point (m„,mrf,ms) = (mo,mo,mo) towards the physical point along 
our m = constant path, the s quark becomes heavier while the u and d quarks 
become lighter. These two effects tend to cancel in any flavour singlet quantity. 
The cancellation is perfect at the symmetric point, and we have found that it 
remains good down to the lightest points we have simulated. 

Since gluonic properties are also flavour singlet, this means that all properties 
of our configurations, from simple ones such as the plaquette, to more complicated 
ones such as the potential and Tq, vary slowly along the trajectory. Compared 
with other paths, the properties of our configurations are already very close to 
those at the physical point. This has many advantages, from technical ones such 
as the rapid equilibration when we move to a new mass point, to physically 
useful results, such as the closeness between partially quenched and full physical 
results. In addition it also enables the lattice spacing to be determined without an 
extrapolation to the physical point, and indeed allows the consistency of various 
definitions to be discussed. 

The flavour symmetry expansion is developed here, by classifying up to 0(5mg) 
how quark-mass polynomials behave under the 5*3 permutation group and the 
5*^7(3) flavour group, leading to Table [3l given for 1 + 1 + 1 quark flavours. 
We also show that for non-chiral (e.g. clover) fermions, where we have different 
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renormalisation for the singlet and non-singlet pieces and also have 0(a) improve- 
ment, that all the additional terms that appear are just these mass polynomials. 
In section 12.51 we classify the hadron mass matrices, and show that certain com- 
binations, for example the Coleman- Glashow relation, have only small violations 
(in terms of the quark mass). 

Turning now to 2 + 1 quark flavours, we have found that the flavour symme- 
try expansion (again when holding the average quark mass, m, constant) leads 
to highly constrained extrapolations (i.e. fits) for non-singlet quantities, such 
as hadronic masses here, and reduce the number of free parameters drastically. 
(There is a short discussion of this point at the end of section 12.31 and in sec- 
tion U31) It is also to be noted that a 2 + 1 simulation is sufficient to determine 
most of the expansion coefficients for the 1 + 1 + 1 case (one exception being the 
particle at the centre of the octet multiplet). 

In section E] we discussed the relationship of the flavour symmetry expansion 
to the chiral perturbation expansion. Lattice simulations are at somewhat large 
pion masses which juxtaposes well with the flavour symmetry expansion presented 
here, while chiral perturbation theory is an expansion about a zero pion mass 
which lattice simulations strive to reach. In section [5] we give, as an example, 
the relationship between these expansions for the pseudoscalar octet. We also 
briefly discuss how a chiral singularity would show up in the flavour symmetric 
expansion and show that at large n the coefficient would drop like a higher power 
in 1/n. (In practice this would be difficult to determine.) 

We have also extended these results in section IHto the partially quenched case 
(when the masses of the valence quarks do not have to be the same as the sea 
quark masses, but we still have the constraint for the sea quarks that fn remains 
constant). In general we have shown (at least to quadratic quark mass order) 
that the number of expansion coefficients does not increase. Thus a (cheaper) 
simulation with partially quenched hadron masses may potentially be of help in 
determining these coefficients. We also show that on the trajectory fn = constant 
the partially quenched error vanishes. 

In sections EHH] numerical results are presented. We first locate in section [9TT] a 
suitable point on the flavour symmetric line, which we take as our initial point for 
the trajectory m = constant. This path can be compared with the trajectory from 
other collaborations. In Fig. [2l] we show the left panel of Fig. [11] again together 
with results from the PACS-CS Collaboration, [IS] and the HS Collaboration, 
[2U] . The strategy of the HS Collaboration was to keep the strange quark mass 
constant (by considering {2M]^-M^)/M^ versus M^/M^, the 'JLAB' plot). We 
see that their points are indeed approximately constant on our plot. 

We then show numerically that flavour singlet quantities, aXs {S = A, A^, 
p and tt) remain constant on the path fn = constant. As the linear term is 
not present, this is a sensitive test of the presence of higher order terms in the 
flavour symmetry expansion and indicates that they appear to be small. This 
also allows an estimation of the scale, a, and a discussion of its consistency. As 
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Figure 24: (2M|, - Ml)/X% versus MljXl for kq = 0.12090. The dashed black 
hne, y = X represents the SU{3) flavour symmetric line. Filled violet circles are on 
32^ X 64 lattices while open violet circles are on a 24^ x 48 sized lattice. Shown are also 
points on the flavour symmetric line (open and filled orange circles). The fits are from 
eq. ([98]) • Results from the PACS-CS Collaboration, [19] and the HS Collaboration, [20] 
are given by cyan coloured squares and magenta coloured diamonds respectively. The 
physical value is denoted by a (red) star. 

can be seen from Fig. [2^ our trajectory does not reach the physical point exactly. 
This is reflected in the fact that different definitions of the scale in Table [15] (last 
column) give results varying by a few percent. 

Results for the hadron mass spectrum are then shown. Numerically we first 
see a mass hierarchy, which confirms our theoretical expectation from the flavour 
symmetric expansion. A series of 'fan' plots for the various multiplets are then 
given, with fits which use the flavour symmetric expansion and show that indeed 
all fits for the pseudoscalar, vector and baryon octets and baryon decuplet are 
highly linear. The higher order terms are very small - one early hint of this is 
the fact that the Gell-Mann-Okubo relations work so well for hadron masses. 
We also note that simulations with a 'light' strange quark mass and heavy 'light' 
quark mass are possible - here the right most points in Figs. [18] - [21] In this 
inverted strange world we would expect the weak interaction decays p — t- S or A. 

In Fig. [25] we plot the results from Table [T7| and compare them with the 
experimental results (also given in this table) for the octet and decuplet hadron 
multiplets. This means that our physical input necessary to determine the hadron 
mass spectrum is kq (i.e. ideally the value corresponding to = m^*), together 
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Figure 25: The masses of the octet and decuplet multiplets as given in Table [T71 using 
Xn to determine the scale, together with the experimental values (short horizontal 
lines) . 

with ml/Xl\* and X^vl*. 

Exploratory partially quenched results for the baryon octet spectrum are 
shown, using the heavier 24^ x 48 data. It is illustrated that they contain useful 
information and allow for the possibility that partially quenched results can help 
in the determination of coefficients in the flavour symmetry expansion. (We plan 
to discuss this further in [T8].) 

We are also applying this method to the computation of matrix elements, [37] , 
some initial numerical results are given in [381 EHl SO] • 
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Appendix 



A The permutation group S3 

If we have three quarks u, d and s with different masses, physics should be un- 
changed if we simply permute the names we give to the quarks. The permutation 
group is not the complete symmetry group - for example we could also perform 
U (1) phase rotations on any particular quark flavour - but it is already enough to 
tell us something useful. The permutation group of 3 objects, 5*3, is the same as 
the symmetry group of an equilateral triangle, Cs^. There are 6 group operations 

1. The identity 

u^u, d ^ d , s — )■ s , (111) 

2. Two cyclic permutations 

u^d, d^s, s^u and u^s, s^d, d^u, (112) 

which correspond to rotations of the triangle through ±120°, and for a 
diagram in the I3-Y plane rotations through ±120°. 

3. Three pair interchanges 

u -H- d , s s; u -H- s , d ^ d; d -(r^ s , u ^ u , (US) 

which correspond to the 3 reflection symmetries of the triangle, and reflec- 
tions of a diagram in the I3-Y plane. 

If an equation is to respect flavour blindness, both sides of the equation should 
transform the same way under all 6 operations. The representations of the group 
allow us to arrange for this to hold. 

The permutation group S3 is a subgroup of SU (3) and has 3 irreducible rep- 
resentations [6]: two different singlets, Ai and A2; and a doublet E. The group 
properties of these are briefly summarised in Table [T] and discussed at greater 
length below. 

A.l Singlet representation Ai 

The representation Ai, the trivial representation, includes objects which are in- 
variant under all 6 group operations. Examples include gluonic quantities, such 
as glueball masses, tq, V{r), as well as certain averages over hadron multiplets. 
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(We shall collectively denote these objects by X.) Examples of quark mass poly- 
nomials with Ai symmetry (complete up to 0(m^)) are 

1 (114) 
m 

, Sm^ + 5171^ + 5ml 

, m((5m^ + + Sml) , SniuSmdSms , 

and linear combinations of these. These are the 7 polynomials of symmetry Ai 
listed in Table [3l the only change is that in the table we have made the re- 
placement m — 7- (m — mo), appropriate for a Taylor expansion about the point 
(mo, mo, TTio). Any other completely symmetric polynomial is a linear combina- 
tion of these, for example 

+ TTT-d + = 3m^ + 3m((5m^ + Sm^ + 5ml) + •^SmuSmdSmg ■ (US) 
A. 2 Singlet representation A2 

This consists of objects which are invariant under cyclic quark permutations 
(triangle rotations), but which change sign under pair exchanges (reflections). 
A2 quantities automatically vanish if any two quark masses are the same. The 
lowest A2 quantity for quark masses is O(m^), 

rfiuml — mdml + m^ms — minis + mlmd — m^mu 
= {6ms - Smu){Sms - Smd){Smu - Smd) . (116) 

Baryon mass combinations with A2 symmetry are 

Mn-Mp- Ms- + Ms+ + Ms- - Mho , (117) 

and the corresponding decuplet quantity, with the p and n replaced by A"*" and 
A°, namely Ma+ - M^o + M^,- - Ms.+ + Mh*o - M^*-. Because particle and 
anti-particle have the same mass, the mesonic analogue of eq. flll7p vanishes. 

Group theory tells us that in a 1 + 1 + 1 flavour world, the splitting, eq. 01171) 
would be proportional to eq. f lll6p and terms of even higher order in m^ (neglect- 
ing electromagnetic effects). 

A. 3 Doublet representation E 

By considering A2 we have found a mass splitting formula for the 1 + 1 + 1 case, 
but by looking at the doublet E we are able to find some more formulae valid for 
the 2 + 1 case which is of more interest for this work. 

The E representation has two states, which mix under the cyclic permutations. 
We can choose to make one state of the doublet even under the reflection u ^ d , 
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and the other state odd. (We could just as well choose any interchange to classify 
our states, but it makes best sense to choose u -H- d, because the hadronic universe 
is almost invariant under that operation.) We have called the even member of 
the doublet E^, the odd member E~ . (There does not appear to be a standard 
notation.) 

An example of an E doublet would be the states 

^ {2\s)-\u)-\d)) and {\u) - \d)) . 



V2 
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It is easily checked that under any group operation they just mix with each other, 
for example under the cyclic operation m — rf, d — )■ s, s — )■ u: 



\u)-\d)) ^ -^{2\u) 

\/3 1 , 



\d) - \s)) 

M))-^^(2|.) 



(119) 



\u) - \d)) 



and so on. In other words, the matrix for a cyclic permutation has the form 



R 



cos 6 =F sin 6 
± sin 6 cos 6 



(120) 



with e = 120°. 

Quark mass terms with E doublet symmetry are 



rrid) 

2\ 



-^(2ms - 

^(2m2 -ml- m^j 



V2 



^(2m3 -ml- ml) 
{limumi + mdml - m^m^ - m„m^) 

-^{muml + m^m^ + 2m^ms - 2m^ms + m^m^ + m„m^) 



1 
1 



rrid) 



^[m^-m^) 
^{mums - mdms 



^{m^ml + mdm 



2m?,m„ 



2m^ms 



mlmd 



m„m. 



:i2ii 



, \{-muml + mdml + mlmd - mum^) } 

The normalisations and phases have been chosen so that each pair transforms 
in the same way as eq. f lllSp under all group operations, i.e. the matrices which 
represent the group operations are the same for every pair. 



B Some group theory 

If the three quarks have equal masses, the QCD Lagrangian is invariant under a 
global U{1) transformation of the quark fields 

V'^e^V, ^^V'e"'^ (122) 
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(corresponding to baryon number conservation) and a global SU (3) flavour trans- 
formation 

ij^Uip, '^^ ^U^ , (123) 

with U a unitary matrix with determinant 1. 

If the quarks are all given different masses we still have the freedom to change 
the phase of each flavour separately, without changing the action, so we have 
conserved currents for each of the three flavours, and three independent U{1) 
symmetries. 

When the quarks have different masses, flavour 5'[/(3) is no longer a symmetry 
of the action, a global SU{2>) rotation no longer leaves the action unchanged, but 
we can still use S'?7(3) to understand the action. 

An analogy of our argument comes from ordinary mechanics or quantum 
mechanics. If we have a quantum mechanical problem which is not rotationally 
symmetric we lose the conservation of angular momentum. But we do still have 
the property that if we rotate the Hamiltonian, H , to give a new problem, with 
the Hamiltonian H' ^ H, then the eigenfunctions of the new Hamiltonian H' can 
be obtained by rotating the eigenfunctions of the original problem. In the case of 
broken flavour symmetry, imposing this (nearly trivial) condition will constrain 
the way in which hadron masses can depend on quark masses. 

Consider the transformation of the quark mass matrix 

M^UMU^ = M', (124) 

(the flavour analogue of a global gauge rotation in colour). The quarks may have 
different masses, Ai' ^ Ai, although they are physically equivalent in the sense 
that the two matrices have the same eigenvalues, but the eigenvectors are rotated 

^' = Uip, ^' = ^U^ . (125) 

Let us now use these deflnitions to investigate the group properties of mass poly- 
nomials. 



B.l Flavour permutations, S3, as a subgroup of SU(3) 

We want to concentrate initially on a set of SU{3) matrices which map a diagonal 
mass matrix to another diagonal matrix when used in eq. ( I124p . These are 

• the identity matrix, 

/ 1 \ 

(126) 
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the cyclic permutations of the quark flavours, 



exp 





exp <^ -i \ -i i \ } (127) 



pair interchanges, 




exp 



exp 




exp < i - 1 1 > (12^ 



Note that when we interchange a quark pair, we also have to change the sign 
of the quarks, to keep the determinant equal to 1, as required for a matrix in 
SU{3). These six matrices are all unitary with determinant 1, so they are all 
members of SU{3). We have also shown that all the matrices can be written 
in the canonical SU{3) form exp{iY^ajXj}. These matrices form a closed set 
under multiplication, with a multiplication table matching that of the group S^,, 
showing that the symmetries of the equilateral triangle are a subgroup of SU{3). 



B.2 Group classification of quark mass polynomials 

This subsection explains how the final column of Table [3] was calculated. 

We can establish many useful results from the S3 subgroup, but it has its 
limitations, it does not connect particles in different permutation sets, see Fig. HI 
By considering 5*3 alone we cannot write down a formula for the mass difference 
between the S° and S~, we cannot even show that the two particles have the 
same mass in the 2 + 1 case. To go further we need to consider the full SU{3) 
group, even though this will involve operations which make the mass matrix 
non-diagonal. 

We can write any SU (3) rotation as a matrix of the form 

= exp |i ^ ajXj I , (129) 
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where the Xj are the 8 Gell-Mann matrices (and aj are real parameters). Here 
we only need to consider infinitesimal transformations 

8 8 

M UMU^ = M + iY,(^j (^j-^ - -^^i) =M+tJ2(^j [^j^M ■ (130) 

i=i i=i 

We write 

Ojip = Xjip 

OjM = [Xj,M] , (131) 

to represent the action of the eight generators of SU{3) on spinors and on matri- 
ces. The eight operators Oj are analogous to the three operators Jj in angular 
momentum. 

In SU (2) we use the eigenvalues of the operator 

J' = i2jl (132) 

i=i 

to identify the irreducible representations of angular momentum. Similarly in 
SU{3) we can use the eigenvalues of the quadratic Casimir operator P [10] 



C = \Y.^I (133) 



to identify irreducible representations of SU{3). (The factor | is a conventional 
normalisation.) Acting on a matrix 

8 8 

C-M = ^ E [^^' t^^' -^]] = i E (^'-^ - + -^^^ • (134) 

We can now begin classifying polynomial functions of Ai. 
At first order, where we have linear functions of mass and Ai can be decom- 
posed as 

8 

M = I ITt[M] + J2 hn^jM] , (135) 

it is simple. We have 

Tt[M] = Mil + M22 + M33 , (136) 
which does not change under SU{3) transformations, so it is singlet. 
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The other elements of Ai can be assigned quantum numbers. For example 
M.21 takes a u quark and changes it to a d, so it has /a = — 1 and hypercharge 
Y = 0. The 6 off-diagonal elements of Ai form the outer ring of the octet, see 
for example Fig. El The two central elements of the octet are the combinations 

2M33-MU-M22 oc Tr[A8A^] 

MU-M22 oc Tr[A3A^]. (137) 

These both have I3 = Y = 0. We can check that both are eigenstates of the 
Casimir operator, with eigenvalue 3, showing that both are pure octet quantities. 
If we make the substitutions 

Mii^niu, M22-^md, M33-^rns, (138) 

we see that the quantities eqs. f ll36l) -( |T371) are proportional to the three linear 
polynomials in Table [3l with the SU{3) assignments given from their behaviour 
when operated on by the Casimir operator. 

It gets more interesting at second order. (Tr[A^])^ and Tr[A^^] are both 
flavour-singlet functions of the mass matrix. It is more convenient to work with 
the linear combinations 

(Tr[A^])' , 3Ti[M^] - iTr[M]f , (139) 

where we have chosen the coefficients so that the second combination will be zero 
at the SU{3) symmetric point. At second order we should be able to construct 
functions of the mass matrix that are in the 1, 8 and 27 representations. One way 
of constructing a quantity that is purely 27-plet is by using the Casimir operator. 
If we take an arbitrary quadratic function of Ai it will usually be a mixture of 
all three representations. If we multiply by 

(C-3)C, (140) 

C will cancel the singlet part, (C— 3) will eliminate the octet part (see Table [TK]) . so 

Representation 1 8 10 TO 27 64 

Casimir eigenvalue 3 6 6 8 15 

Table 18: The eigenvalues of the quadratic Casimir operator, C, eq. (jl33p . for the 
SU (3) representations needed in this article. 

the operator eq. f ll40p leaves a pure 27-plet function of Ai. Using the eigenvalues 
in Table [18] we can construct similar operators to project out objects in the other 
representations of SU{3). Of course it would be tedious to do this by hand: we 
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have programmed the group operations in Mathematica so that the group theory 
can be done more easily and rapidly. 

Another useful technique is to use the raising and lowering operators I±, U±, 
V± [l2] to move around within a multiplet. Once we have one state in a multiplet, 
these operators allow us to construct all the other states. Because infinitesimal 
SU{3) operations do not preserve diagonality, a typical eigenstate of the Casimir 
operator will involve all nine elements of the quark mass matrix Ai, not just the 
three diagonal elements. For example, if we explicitly write out the SU{3) singlet 
quantity 3Tr[A^2] _ (Tr[A^])2 in eq. ( 1T391) it is 

Pi = 2MnMii + 2M22M22 + 2M33M33 

+ 6A^i3A^3i + 6M23M32 (141) 

-2MuM22 - 2A^iiA^33 - 2A^22-M33 • 

We can use the techniques discussed in this section to write down a pure SU{3) 
27-plet quantity, with the same 5*3 properties as eq. f ll4ip : the result is 

Pat = MuMu + M22M22 + M33M33 

-M12M21 - M13M31 - M23M32 (142) 
-M11M22 - M11M33 - M22M33 . 

Expressed this way, the 27-plet and singlet are clearly different functions of the 
full 9-element mass matrix. However, if we just consider a diagonal mass matrix, 
Aiij = if i 7^ j, Aiu = niu, M.22 = "^d, A^33 = then the quantities become 
indistinguishable: 

Pi 2{ml + ml + ml - rriumd - mums - mdmg) 
= 3(5m^ + Sm"^ + 5ml) 
P^^ ml + ml + ml - muma - mums - mams 

= ^(5ml + 6mUSml). (143) 

Both collapse to the same quark mass polynomial, Sm^ + Sm'^ + Sm^, so this 
polynomial is allowed to appear in equations for singlet and 27-plet physical 
quantities, but not in equations for any other SU{3) representation. This poly- 
nomial is recorded in Table [3] with the symmetry representations Ai and 1 or 
27. 

We can continue and use the methods of this subsection to classify all poly- 
nomials up to cubic order, the results are recorded in Table O 

B.3 Matrix representations of SU(3) 

To construct hadron mass matrices for octet and decuplet hadrons we need to 
analyse 8x8 and 10 x 10 matrices by their 5*3 and SU{3) properties. 
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To get started we need to construct 8x8 and 10 x 10 representations of the 
SU{3) generators. We can do this by considering the known behaviour of the 
hadron multiplets under the SU (2) subgroups, isospin J, [/-spin and \^-spin, and 
the hypercharge, Y, [i2]: 

Ai = 2Ji A2 = 2/2 A3 = 2/3 

A4 = 2Vi As = 2V2 (144) 
Ae = 2Ui A7 = 2f/2 As = 73 y . 

These 8x8orl0xl0A matrices have the same commutation relations as the 
usual 3x3 matrices 

[A„ Xj] = 2ifijk\u . (145) 

Once we have the eight A matrices for our hadron multiplet we can use eq. (I134p 
to classify any other matrices. For the hadron mass matrices, we need all the 
flavour-conserving matrices. For the decuplet mass matrix these are all diago- 
nal matrices; but for the octet mass matrix they can include some off-diagonal 
elements, because the E° and A have the same flavour quantum numbers. Our 
results for the decuplet and octet matrices are given in Tables H] and [51 

We have other methods of constructing the matrix representations of S'?7(3). 
In addition to the Casimir projection method sketched here, we can start with 
one matrix which belongs to a known SU{?>) representation, and then build all 
the other matrices in that representation by repeatedly acting with raising and 
lowering operators. For example, in the decuplet case we know that the 10 x 10 
matrix 
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must be a pure 64-plet, because it interchanges the A~ and the which changes 
strangeness by ±3. From Fig. [7] we see that the 64-plet is the only representation 
in 10 ® 10 that can change strangeness by 3 units. Starting from the matrix 
in eq. fll46p we can construct a set of 64 matrices which transform amongst 
themselves under all group operations eq. fll3ip . 

Once we have classified all the 10 x 10 and 8x8 matrices according to their 
SU (3) and 5*3 behaviour, we can read off the rows of Tables [Hand [51 For example, 
knowing that the following 8x8 matrix is an octet with symmetry E~ gives the 
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fifth row of Table [5l 



/ -1 



000 000 0\ 












1 



75 

0^ 



1 



(147) 



V 



00 00-1/ 



B.4 Hadron mass matrices 

We describe the hadron masses via a hadron mass matrix "H, a 10 x 10 matrix 
for the decuplet baryons, an 8 x 8 matrix for octet baryons or mesons. 

If the different quark flavours have different masses, a global SU{3) rotation 
of the quark mass matrix, eq. fll24p . leads to a change in the quark mass matrix, 
Ai — )• A^', but does not change the eigenvalues of the matrix, or the essential 
physics of the situation. 

What will be the effect of a flavour rotation of the quark Lagrangian on a 
hadronic mass matrix 

If we are considering a Taylor expansion for hadronic masses, all the elements 
of Ti will be polynomials of the elements in the quark mass matrix, Aiij. Flavour 
blindness requires that we still get equivalent physics when we change J^ij — )■ 
■M'ij, i.e. that the eigenvalues of Ti are unchanged, and the eigenvectors of Ti 
rotate according to eq. (11251) . Writing this as an equation. 



Using the unitarity of U we can rewrite this as an invariance condition. 



The effect of changing Ai to M.' can be exactly cancelled by the effect of an 
SU{3) rotation on Ti. 

To construct an invariant matrix satisfying eq. (I149p we have to pair up ma- 
trices of known symmetry, constructed as described in Appendix IB.3t with A4 
polynomials of known symmetry, constructed using the methods of Appendix lB.2[ 
This gives us a hadron mass matrix of the form 



n' = n{M[A = UH{M^j)U^ . 



(148) 



(149) 



n 




(150) 
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To give an SU{2) analogy, we can make a rotationally invariant system (i.e. a 
system with total spin zero), by coupling together two particles with the same 
J, but not by coupling together two particles with different J. Similarly, to give 
a hadron mass matrix under the SU{3) operation, eq. f ll49p . we must give every 
matrix a coefficient of the same symmetry, as shown schematically in eq. (11501) . 

Once we have (with the help of Mathematica) , constructed the most general 
matrix satisfying eq. (11491) . we make the substitutions A^n — )■ A^22 — > ^d, 
Ai33 ^ rUs, and A^jj — )■ if i 7^ j to get mass formulae for all the hadrons. 

We now consider an example. In TableOwe list the 6 matrices which can occur 
in the octet meson mass matrix in the 1 + 1 + 1 flavour case. In the 2 + 1 flavour 
case (m„ = m^) the two E" matrices drop out, because their coefficients must be 
odd under the exchange u -H- d, leaving just 4 matrices which can contribute. We 
calculate the most general form of the meson mass matrix, by demanding that it 
is invariant under SU{3) rotations, eq. (I149p . and flnd that in the 2 + 1 case with 
m = constant we get 
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keeping terms up to quadratic order. From this we read off 



Ml 
Ml 
Ml 



Ml - 2as5mi + (hi - 2^8 
Ml + asSmi + {h + + 
Ml + 2a85mi + (6i + 2^8 



— 3b27)6mi 



27b27)5mt 

We can checls: that these equations are consistent with Table 



3Ml + AMI + Ml 
-3Ml + 2Ml + Ml 
-Ml + AMI - 3Ml 



8{Ml + bi5m1) 
10(a85m; + bs5mf] 
120b27S'ml . 



(151) 



(152) 



(153) 



An alternative method, as discussed in section [31 would be to start from the 
simultaneous equations in eq. fll53p and solve the system to derive eq. fll52p . 
Note that in eq. ( H3|l we have also re-written the results in a form to agree with 
the notation of the partially quenched results, so as = —a and 



h 

bs 

b27 



/3o + 4/3i + 6/32 + 1/33 
/3i + 3/32 + 



-^/33 



(154) 



C Coordinate choice for partially quenched for- 
mulae 

It is often convenient to plot quantities against the pseudoscalar meson mass 
squared, because then we know better the location of the physical point and the 
chiral limit. If we do want to use pseudoscalar mesons, the best choice is to 
replace the light sea quark by the full (non partially quenched) pion, the light 
valence quark by the partially quenched pion, and to replace the valence strange 
quark mass by the partially quenched s^aiSvai meson (the 'strange pion'), which 
we call the rjs. This is a particle that doesn't exist in the real world, but which 
we can measure in the partially quenched channel. Determining the valence s 
quark mass from the kaon has disadvantages, as we shall shortly see. 



76 



We introduce the mesonic variables 



X = Mlf^a -M% = 2a5mi + l3o5mf + 2(3i5m^ 

z = Ml-M% = 2a5^s + PoSmf + 2Pi6fil, (155) 

keeping terms up to second order in the quark masses. In terms of these variables 
the decuplet mass formulae eq. ( l53l) become 

Ma = Mo + 3Ay + Box^ + 3Biy^ 

Ms* = Mo + A{2y + z) + Box^ + B^{2y^ + z^) + B^iz ~ yf 

Ms* = Mo + A{y + 2z) + Box^ + Bi{y^ + 2z^) + B2{z ~ 

Mn = Mo + 3Az + Box^ + 3Biz^ , (156) 



with 



A 

Bo 
Bi 



A 
2^ 

2aBo - 3Al3o 



aBi - A/3i 



The form of eq. f ll56p exactly repeats the form of eq. ( l53l) . but the new constants 
involve a combination of curvature terms from the pion mass equation and from 
the baryon mass equation. 

Suppose we use the PQ kaon mass (instead of the strange pion) to represent 
the strange quark mass, i.e. we replace z defined in eq. fll55p by 

w = 2M|^pg - M^PQ - Mli^u 

= 2a5fXs + l3o5mf + 2l3i5fxl + 2(32{5fXs- Sfxif . (158) 

At first order, w is just as good as z, but if we are interested in curvature, it is less 
suitable, because at second order it involves both the valence s and the valence 
/, unlike eq. f ll55p . Using w instead of z, the decuplet mass formulae become 

Ma = Mo + 3Ay + I3ox^ + 3Biy^ 

Ms, = Mo + A{2y + w) + Box^ + B,{2y^ + w^) + B2{w - yf + Bxiw - 

Ms. = Mo + A{y + 2w) + Box^ + Bi{y^ + 2w^) + B2{w-y)^ + 2Bx{w-yy 

Mn = Mo + 3Aw + BoX^ + 3BiW^ + 3Bxiw-yY, (159) 

with A, Bo, Bi, B2 defined as in eq. fll57p . but with an extra curvature coefficient 

^x = -§. (160) 
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so one fit constraint is lost (or deeply hidden) if we use the kaon mass to represent 
the strange mass. 

Finally, we want to relate the partially quenched fit to the unitary results, on 
our trajectory ^{2mi + ms) = mo. If we use bare quark masses as our coordinates, 
we do this by using the substitutions 

6fii — 7- 6mi , 6fis — !■ —26mi , (161) 

giving 

Ma = Mo + 3A5mi + [Bq + 3Bi]6m^ 
Ms* = Mo + [Bo + 6Bi + 9^2] 5mf 

Ms* = Mo- 3A6mi + [Bq + 9Bi + 9B2]Sm^ 

Mn = Mo- 6A6mi + [Bo + 12Bi]6m^ . (162) 

However, if we use meson-based coordinates, such as eq. f llSSp . the mapping back 
to the unitary result is more complicated. 



y — )■ X 

z — )■ — 2x H X . 163 

The mapping from 2;, our measure of the strange quark mass, back to x is com- 
plicated by a second order term. The reason is clear. On our trajectory, the 
relation 2bmi + = is made exactly true for bare lattice quark masses, while 
the meson mass relations 2M^ + M'i fa constant or 2M?- + M^ constant are 
only true to leading order. Thus in conclusion if we are considering the curvature 
terms it is definitely better to use (bare) lattice quark masses as the coordinates. 



D The action 

The particular clover action used here has a single iterated mild stout smearing, 
j43] for the hopping terms together with thin links for the clover term (this ensures 
that the fermion matrix does not become too extended). Together with the (tree 
level) Symanzik improved gluon action this gives 

S = Spu + Spd + Sps , (164) 

with the gluon action 




(165) 
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and 

R 6co 10 20 1 

P = -?r = ^ and Co = — , ci = -— . (166) 
i/O i/o -"-^ -"-^ 

For each flavour the Wilson-Dirac fermion action is 



S 



Fq 




\ X][?(a;)(7M - l)Ut,{x)q{x + afx) - g(x)(7^ + l)t/J(x - a/i)g(x - a/i)] 
— g(a;)g(x) - ^ac^^ ^ g(x)o-^j.F^^(a;)g(x) I , (167) 



2 



2Ka ' 4 



where F is the 'clover' field strength, necessary for 0(a)-improvement. As the up 
and down quarks are always taken here as mass degenerate we have Ku = k^, = ki- 
To keep the action highly local, the hopping terms use a stout smeared link 
('fat link') with a = 0.1 'mild smearing' for the Dirac kinetic term and Wilson 
mass term, 



= I [V,Ul - U,Vl - \Tv{V,Ul - U.Vi)] , (168) 

where V^(a;) is the sum of all staples around f/^(x). The clover term is built from 
thin links as it is already of length 4a and, as previously mentioned, we do not 
want the fermion matrix to become too extended. Stout smearing is analytic and 
so a derivative can be taken (so the HMC force is well defined) and also allows 
for perturbative expansions |44j . 

The clover coefficient, c^u,, has recently been non-perturbatively fixed, [21], 
by requiring that the axial Ward identity (WI) quark mass determined in several 
different ways is the same. A sensitive way of achieving this is the Schrodinger 
functional formalism. Further details of our results may be found in [23]. Csw 
is determined for 3 mass degenerate or SU{'i) flavour symmetric quarks (where 
Ki = Kg = kq) in the chiral limit. A 5th order polynomial in interpolating 
between the numerically determined Cswido) points was found to be 



cl^{go) = 1 + 0.269041 + 0.29910 g'^ - 0.11491 g^ - 0.20003 g^ + 0.15359 g^^ . 

(169) 

(This interpolation function is constrained to reproduce the 0((7q) perturbative 
results, [44j, in the /3 — )■ oo limit and therefore has four free fit parameters.) We 
take this result to define Csw for a given (3. 

Improving one on-shell quantity to O(a^) (here the axial WI quark mass) fixes 
Cswigi) then all masses are automatically improved to O(a^), 

Mm Mm 
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rather than just to 0{a). Operators in general require further 0(a) operators 
together with associated improvement coefficients to ensure 0(a)-improvement 
for physical on-shell quantities. 

This determination of Csw via the Schrodinger functional formalism also pro- 
vides an estimate for the critical kq, |2l], of 

«0;c(^o) = \ [l + 0.002391 c/^ + 0.0122470 - 0.0525676 
8 

+0.0668197^^- 0.0242800^0^°] . (171) 

(Again this interpolation function is constrained to reproduce the 0((7q) pertur- 
bative results, [H], in the /3 — )■ oo limit. The errors for c*^ from the fit are 
estimated to be about 0.4% while for re* we have 0.02% at /3 = 14.0 rising to 
0.15% at /3 = 5.10.) 

The simulations only need knowledge of Csw to proceed; however it is useful to 
check consistency between different determinations of fi;o;c (via the Schrodinger 
functional or the pseudoscalar mass). For /3 = 5.50 then using eq. f ll7ip we find 
Ko;c = 0.120996 (the direct simulation result is kq.^ = 0.121125(330), |2l]). This 
is to be compared with the estimation in section 18.21 which is quite close. (It 
should also be noted that different determinations should only agree up to O(a^) 
effects.) 

E Hadron masses 

We collect here in Tables [T9] - [271 values of the pseudoscalar octet, vector octet, 
baryon octet and baryon decuplet masses. In Table [19] we give values along the 
flavour symmetric line {ki = Kg = kq), while in Tables [20] - [23] and in Tables 12^ - 
[27] we give results for kq = 0.12090 and kq = 0.12092 respectively, while keeping 
m = constant, eq. (IHSj) . 

In Tables [28] - [30] we give the ratios (i.e. hadron octet or decuplet masses 
normalised with their centre of mass). 

The data sets are roughly ~ 0(2000) trajectories for the 24^ x 48 lattices and 
0(1500) - 0(2000) trajectories for the 32^ x 64 lattices (with the exception for 
the Kq = 0.12095 results which are ~ 0(500) trajectories). The errors are all 
taken from a bootstrap analysis of the ratio (which often enables a smaller error 
to be given for the ratios than simply using error propagation). 
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Ko 


Nl 


X 


Nt 




aMp aMN 




0, 


,12000 


16^ 


X 


32 


0.4908(17) 


0.6427(23) 0.9612(42) 


1.048(6) 


0, 


,12030 


16^ 


X 


32 


0.4026(19) 


0.5635(38) 0.8374(74) 


0.9414(107) 


0, 


,12050 


243 


X 


48 


0.3375(24) 


0.4953(47) 0.7201(83) 


0.8216(89) 


0. 


,12080 


243 


X 


48 


0.2260(10) 


0.3903(55) 0.5417(68) 


0.6415(99) 


0. 


,12090 


16^ 


X 


32 


0.2209(49) 


0.4192(97) 0.6298(251) 


0.7811(274) 


0, 


,12090 


243 


X 


48 




see Tables [20] -[21 




0, 


,12090 


32^ 


X 


64 




see Tables [20] -[21 




0, 


,12092 


243 


X 


48 




see Tables [21 -[23 




0, 


,12095 


32^ 


X 


64 


0.1508(4) 


0.3209(27) 0.4329(49) 


0.5541(80) 


0, 


,12099 


32^ 


X 


64 


0.1297(10) 


0.3154(67) 0.4127(117) 


0.5476(168) 



Table 19: The results for the hadrons on the symmetric line, aM^^, aMp, aM^ and 
oMa for (/3,Cs„,a) = (5.50,2.65,0.1). 



K 



aM„ 



16-' X 32 



(0.121040, 0.120620) 0.1962(74) 0.2447(49) 0.2773(37) 



24^ X 48 



(0.120830, 0.121040) 
(0.120900, 0.120900) 
(0.120950, 0.120800) 
(0.121000, 0.120700) 
(0.121040, 0.120620) 



0.1933(6) 

0.1779(6) 

0.1661(8) 

0.1515(10) 

0.1406(8) 



0.1688(7) 
0.1779(6) 
0.1845(7) 
0.1898(8) 
0.1949(6) 



0.1391(11) 

0.1779(6) 

0.2011(7) 

0.2209(6) 

0.2361(5) 



32^ X 64 


(0.120900, 0.120900) 
(0.121040, 0.120620) 
(0.121095, 0.120512) 
(0.121145, 0.120413) 


0.1747(5) 
0.1349(5) 
0.1162(8) 
0.09694(88) 


0.1747(5) 
0.1897(4) 
0.1956(5) 
0.2016(4) 


0.1747(5) 
0.2321(3) 
0.2512(3) 
0.2683(3) 



Table 20: The results for the pseudoscalar octet mesons: aM^, aMx and aM^^ 
{|3,Cs^,,a) = (5.50,2.65,0.1) where = 0.12090. 



for 
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{ki,Ks) aMp aMx* aM^^ 

16^ X 32 



(0.121040, 0.120620) 0.4353(123) 0.4331(84) 0.4380(60) 



24^ X 48 


(0.120830, 0.121040) 
(0.120900, 0.120900) 
(0.120950, 0.120800) 
(0.121000, 0.120700) 
(0.121040, 0.120620) 


0.3460(22) 
0.3494(25) 
0.3400(40) 
0.3364(43) 
0.3270(50) 


0.3335(30) 
0.3494(25) 
0.3473(32) 
0.3517(30) 
0.3484(28) 


0.3198(48) 
0.3494(25) 
0.3546(27) 
0.3663(20) 
0.3701(18) 


32^ X 64 


(0.120900, 0.120900) 
(0.121040, 0.120620) 
(0.121095, 0.120512) 
(0.121145, 0.120413) 


0.3341(34) 
0.3127(38) 
0.3123(43) 
0.3210(63) 


0.3341(34) 
0.3380(21) 
0.3426(20) 
0.3500(24) 


0.3341(34) 
0.3632(14) 
0.3738(11) 
0.3880(11) 



Table 21: The results for the vector octet mesons: aMp, aMx* and aM<^^ for 
(/?, Csyj,a) = (5.50, 2.65, 0.1) where kq = 0.12090. 



{ki,Ks) oMjv aMA aMs aM^ 

16^ X 32 



(0.121040, 0.120620) 


0.5817(214) 0.5941(182) 


0.6311(128) 


0.6353(121) 


24=^ X 48 


(0.120830, 0.121040) 

(0.120900, 0.120900) 
(0.120950, 0.120800) 
(0.121000, 0.120700) 
(0.121040, 0.120620) 


0.4976(25) 

0.4811(33) 
0.4737(68) 
0.4648(46) 
0.4466(66) 


0.4859(43) 

0.4811(33) 
0.4794(58) 
0.4815(49) 
0.4810(57) 


0.4791(31) 

0.4811(33) 
0.4871(55) 
0.4910(36) 
0.4843(42) 


0.4679(39) 

0.4811(33) 
0.4938(48) 
0.5055(28) 
0.5068(32) 


32=^ X 64 


(0.120900, 0.120900) 
(0.121040, 0.120620) 
(0.121095, 0.120512) 
(0.121145, 0.120413) 


0.4673(27) 
0.4267(50) 
0.4140(61) 
0.4016(89) 


0.4673(27) 
0.4547(43) 
0.4510(58) 
0.4507(65) 


0.4673(27) 
0.4697(33) 
0.4690(37) 
0.4761(39) 


0.4673(27) 
0.4907(21) 
0.4971(21) 
0.5092(19) 



Table 22: The results for the octet baryons: aM^-, aM\, aM-^ and aMs for 
{I3,csw,a) = (5.50,2.65,0.1) where kq = 0.12090. 
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aMa 




16 


3 X 32 






(0.121040, 0.120620) 


0.7437(227) 


0.7490(184) 


0.7537(146) 


0.7595(114) 


24^ X 48 


(0.120830, 0.121040) 


0.5906(73) 


0.5801(89) 


0.5685(114) 


0.5548(151) 


(0.120900, 0.120900) 


0.5933(88) 


0.5933(88) 


0.5933(88) 


0.5933(88) 


(0.120950, 0.120800) 


0.5817(55) 


0.5895(48) 


0.5973(43) 


0.6050(38) 


(0.121000, 0.120700) 


0.5883(101) 


0.6006(77) 


0.6133(61) 


0.6262(51) 


(0.121040, 0.120620) 


0.5483(137) 


0.5679(90) 


0.5902(64) 


0.6108(48) 




32 


^ X 64 






(0.120900, 0.120900) 


0.5675(64)) 


0.5675(64) 


0.5675(64) 


0.5675(64) 


(0.121040, 0.120620) 


0.5520(79) 


0.5744(48) 


0.5968(34) 


0.6194(28) 


(0.121095, 0.120512) 


0.5161(185) 


0.5541(98) 


0.5812(52) 


0.6104(33) 


(0.121145, 0.120413) 


0.5071(211) 


0.5576(105) 


0.6018(51) 


0.6420(29) 



Table 23: The results for the decuplet baryons: aM^, aM^*, aM^* and aMfi for 
{P,Csw,a) = (5.50,2.65,0.1) where kq = 0.12090. 



{ki, Kg) aM„ aMx aM^^ 

24=^ X 48 

(0.120920, 0.120920) 0.1694(9) 0.1694(9) 0.1694(9) 

32^ X 64 

(0.121050, 0.120661) 0.1280(6) 0.1813(5) 0.2221(4) 



Table 24: The results for the pseudoscalar octet mesons: aM^^, aMx and aiW^^ for 
{(3,Csw,C() = (5.50,2.65,0.1) where kq = 0.12092. 



{ki, Ks) aMp gMk* aM^^ 

24^ X 48 

(0.120920, 0.120920) 0.3404(44) 0.3404(44) 0.3404(44) 

32^^ X 64 

(0.121050, 0.120661) 0.3161(38) 0.3354(22) 0.3564(16) 



Table 25: The results for the vector octet mesons: aMp, qMk* and aM^^ for 
{f5,Csnj,Oi) = (5.50,2.65,0.1) where kq = 0.12092. 
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{ki,Ks) uMn oMa aMs aMw 





243 


X 48 






(0.120920, 0.120920) 


0.4725(39) 


0.4725(39) 


0.4725(39) 


0.4725(39) 




323 


X 64 






(0.121050, 0.120661) 


0.4127(42) 


0.4444(35) 


0.4580(31) 


0.4798(22) 



Table 26: The results for the octet baryons: aMj^, aM\, aM^, and aMs for 
{P,Csy,,a) = (5.50,2.65,0.1) where = 0.12092. 



(«;,«s) 






aMs* 


aMn 




243 


X 48 






(0.120920, 0.120920) 


0.5790(97) 


0.5790(97) 


0.5790(97) 


0.5790(97) 




32^ 


X 64 






(0.121050, 0.120661) 


0.5457(108) 


0.5607(72) 


0.5800(51) 


0.6005(40) 



Table 27: The results for the decuplet baryons: oMa, aMs*, aMs* and aMo, for 
{P,Csy„a) = (5.50,2.65,0.1) where = 0.12092. 











243 X 48 


(0.120830, 0.121040) 
(0.120900, 0.120900) 

(0.120950, 0.120800) 
(0.121000, 0.120700) 
(0.121040, 0.120620) 


1.025(2) 
1.0 

0.9859(22) 
0.9706(34) 
0.9581(60) 


0.9877(12) 
1.0 

1.007(1) 
1.015(2) 
1.021(3) 


0.9470(155) 
1.0 

1.028(6) 

1.057(12) 

1.086(12) 


32^ X 64 


(0.120900, 0.120900) 

(0.121040, 0.120620) 
(0.121095, 0.120512) 
(0.121145, 0.120413) 


1.0 

0.9488(50) 
0.9392(63) 
0.9431(109) 


1.0 

1.026(3) 
1.030(3) 
1.028(5) 


1.0 

1.102(6) 
1.124(7) 
1.140(9) 



Table 28: Ratio results for the vector octet mesons: Mp/Xp, MK*/Xp and M^^/Xp 
for (/?, Csw,oi) = (5.50, 2.65, 0.1) where kq = 0.12090. 
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Mn/Xn 




M^/Xn 


Ms/Xn 


X 48 


(0.120830, 0.121040) 


1.033(2) 


1.009(6) 


0.9949(13) 


0.9717(26) 


(0.120900, 0.120900) 


1.0 


1.0 


1.0 


1.0 


(0.120950, 0.120800) 


0.9769(33) 


0.9887(84) 


1.005(1) 


1.018(3) 


(0.121000, 0.120700) 


0.9543(32) 


0.9885(77) 


1.008(2) 


1.038(3) 


(0.121040, 0.120620) 


0.9319(56) 


1.004(7) 


1.011(2) 


1.058(4) 




32^ 


X 64 






(0.120900, 0.120900) 


1.0 


1.0 


1.0 


1.0 


(0.121040, 0.120620) 


0.9229(47) 


0.9833(58) 


1.016(2) 


1.061(4) 


(0.121095, 0.120512) 


0.8999(77) 


0.9804(111) 


1.019(4) 


1.081(8) 


(0.121115. 0.120113) 


0.8688(118) 


0.9919(130) 


1.030(5) 


1.101(8) 



Table 29: Ratio results for the octet baryons: Mj^/Xj^, M\/X]^, My,/Xj^ and 
M^/Xn for (/3,c,^,a) = (5.50,2.65,0.1) where acq = 0.12090. 







Me./Xa 


MsV^A 


Mn/^A 




243 


X 48 






(0.120830, 0.121040) 


1.021(6) 


1.003(2) 


0.9824(44) 


0.9588(121) 


(0.120900, 0.120900) 


1.0 


1.0 


1.0 


1.0 


(0.120950, 0.120800) 


0.9868(14) 


1.0002(3) 


1.013(2) 


1.026(3) 


(0.121000, 0.120700) 


0.9790(42) 


0.9993(24) 


1.020(6) 


1.042(8) 


(0.121040, 0.120620) 


0.9634(72) 


0.9978(53) 


1.037(11) 


1.073(14) 




323 


X 64 






(0.120900, 0.120900) 


1.0 


1.0 


1.0 


1.0 


(0.121040, 0.120620) 


0.9609(44) 


0.9999(28) 


1.039(7) 


1.078(9) 


(0.121095, 0.120512) 


0.9426(120) 


1.012(10) 


1.062(19) 


1.115(24) 


(0.121145, 0.120413) 


0.9185(145) 


1.010(13) 


1.090(23) 


1.163(29) 



Table 30: Ratio results for the decuplet baryons: Ma/X^, Ms*/Xa, Ms*/Xa and 
Mu/Xa for {l3,Csyj,a) = (5.50,2.65,0.1) where kq = 0.12090. 
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